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Honey  bees  (Apis  mellifera  L.)  of  African  and  European  ancestry  were 

distinguished  by  analysis  of  restriction  fragment  length  polymorphisms  (RFLPs) 

defined  by  two  enzymes  at  a  highly  polymorphic  locus  corresponding  to 

genomic  probe  pB178.      Thirty-six  Msp\   restriction  fragment  patterns,   or 

variants,  and  thirty-three  Dde\  variants  were  identified.     Sixty-five  pairwise 

combinations  of  the  Msp\  and  Dde\  variants,  referred  to  as  alleles,  were  found 

among  the  individual  drones  tested.   Variants  and  alleles  were  discontinuously 

distributed  in  USA  and  South  African  drones:    only  one  Dde\  variant  and  no 

Msp\  variant  (hence  no  allele)  were  common  to  these  two  sample  populations. 

I 
The  diversity  in  the  variants  and  alleles  found  in  the  South  African  drones  was 
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greater  than  in  the  USA  drones.  Msp\  variants  were  discontinuously  distributed 
among  workers  bees  of  the  subspecies -4.  m.  mellifera  (west  European),  A.  m. 
ligustica  and  A.  m.  caucasica  (east  European),  and  A.  m.  scutellata  (South 
African). 

Ancestry  in  New  World  bees  was  inferred  from  variant  and  allele 
frequencies  at  this  locus.  In  USA  bees,  variants  characteristic  of  east  and  west 
European  bees  were  found  at  frequencies  consistent  with  previously  identified 
nuclear  and  mitochondrial  DNA  markers.  In  neotropical  drones,  European  and 
African  origins  were  assumed  for  variants  that  were  similar  in  fragment  pattern 
to  variants  in  the  USA  and  South  Africa,  respectively.  The  alleles  identified  in 
the  majority  of  neotropical  colonies  were  African.  There  was  little  evidence  of 
hybridization  of  African  and  east  European  bees  in  neotropical  bees:  variants 
specific  to  South  African  bees  were  detected  at  the  highest  frequencies,  while 
east  European-specific  variants  were  absent  or  detected  at  very  low 
frequencies.  A  variant  that  may  be  specific  for  A.  m.  mellifera  vjas  found 
throughout  the  neotropics  at  frequencies  that  varied  from  1 9%  to  33%  but  did 
not  increase  clinally  in  a  northward  direction.  ;  v^ 

Determination  of  the  allelic  relationships  between  Msp\  fragments  and 
between  Dde\  fragments  was  initiated  by  mapping  these  restriction  sites  within 
the  probe,  and  correlating  fragments  on  Southern  blots  to  the  locations  of 
restriction  sites  on  the  resulting  map.  "     ' 
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INTRODUCTION 

Descendants  of  ancestral  honey  bees  migrated  from  the  site  of  their 

■J  ■■' 

origin,  presumed  to  be  northeast  Africa  or  the  Middle  East,  and  became  widely 

distributed  in  the  Old  World  (Garnery,  Cornuet  &  Solignac  1992).   Subsequent 

adaptation  to  a  wide  variety  of  ecological  conditions  resulted  in  the  evolution 

of  subspecies  of  the  honey  bee.  Apis  mellifera  L.   (Ruttner  1988).     These 

subspecies  are  distinguished  on  the  basis  of  physical,  behavioral,  and  ecological 

criteria  and  are  defined  quantitatively  by  discriminant  analysis  of  morphological 

characters  (Daly   &  Balling    1978;   Ruttner   1988;    Ruttner,   Tassencourt  & 

Louveaux  1978). 

Honey  bees  were  first  introduced  to  the  New  World  by  17th  and  18th 

century  European  settlers.    For  over  two  hundred  years,  the  genetic  diversity 

found  in  New  World  honey  bees  resulted  from  the  importation  of  primarily  four 

European  subspecies  or  races:     west  and  north  European  A.  m.  mellifera 

(referred  to  here  as  west  European);  north  Mediterranean y4.  m.  ligustica  Spinola 

and  A.  m.  carnica  Pollmann  (which  in  this  report  will  be  called  east  European; 

Ruttner  1 988);  and  the  east  European  bee  A.  m.  caucasica  Gorbachev  (referred 

to  here  as  east  European)  (Kent  1988;  Kerr,  DeLeon  &  Dardo  1982). 


-nK  •'■~'  —  -'~  >■'•'-  ■■ —     ■■ —    -'' -'   - 


A," 


.V>' 


Despite  abundant  nectar  sources  in  the  neotropics,  self-sustaining  feral 

populations  of  European  bees  did  not  become  established  (Michener  1975). 

The  poor  performance  of  European  bees  in  the  tropics  has  been  attributed  to 

I 
their  failure  to  adapt  to  environmental  conditions  vastly  different  from  those  in 

which  they  evolved  (Ruttner  1 988).  To  improve  commercial  honey  production, 

queens  of  a  central  and  south  African  race,  A.  m.  scutellata  Lepeletier,  were 

imported  into  Brazil  in  1956  for  experimental  breeding  with  European  bees. 

Swarms  of  African  bees  escaped  into  the  wild  (Kerr  1967)  and  subsequently 

proliferated  to  form  large  feral  populations  that  spread  through  tropical  South 

and  Central  America.  African  bees  entered  Texas  in  1 990  and  Arizona  in  1 993. 

I  . 

The  release  and  spread  of  African  bees  has  been  disruptive  to  the 

behavior  and  ecology  of  the  melliferous  flora  and  fauna  in  the  neotropics 

(Cantwell  1974;  Michener  1975;  Roubik  1980,  1989;  Spivak,  Fletcher  &  Breed 

1991;  Taylor  1977;  Winston,  Taylor  &  Otis  1983).     Explanations  for  the 

migration  pressure  resulting  from  the  introduction  of  the  bees  from  Africa  and 

for  the  rapid  colonization  of  New  World  territories  previously  unoccupied  by 

honey  bees  have  focused  on  the  extent  to  which  neotropical  bees  represent  an 

:/,  I 

admixture  of  African  and  European  subspecies  (Hall  1990,  1991,  1992a;  Hall 

&  Muralidharan  1989;  Rinderer  1986;  Rinderer  ef  a/.  1985,  1991;  Sheppard  ef 

al.  1 991 ;  Smith,  Taylor  &  Brown  1 989;  Taylor  1 985).  The  ability  to  distinguish 

honey  bee  subspecies  is  essential  to  ascertain  the  extent  of  hybridization.    In 

turn,  the  degree  of  hybridization  will  influence  the  northern  limit  of  African  bee 
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introgression  in  the  USA,  about  which  there  is  considerable  debate  (Dietz  1 986; 
Dietz,  Krell  &  Pettis  1986;  Rinderer  1986;  Roublk  1986;  Southwick,  Roubik  & 
Williams  1990;  Taylor  1977;  Taylor  &  Spivak  1984;  Villa,  Gentry  &  Taylor 
1987;  Villa,  Rinderer  &  Collins  1993). 

The  retention  of  African  morphology  and  behavior  in  neotropical  bees  has 
stimulated  the  quest  for  the  identification  of  their  ancestry,  which  will  aid  in 
understanding  the  mechanisms  contributing  to  their  phenomenal  success. 
Efforts  to  ascertain  the  extent  of  interaction  between  the  extant  European  bees 
and  the  African  bees  in  the  neotropics  have  been  hindered  by  a  lack  of  genetic 
markers  specific  for  each  Old  World  subspecies  introduced.  Intermediate  values 
of  morphometric  values  and  allozyme  frequencies,  taxonomic  characters 
traditionally  used  to  identify  honey  bee  subspecies  and  subspecies  groups,  can 
only  suggest  that  hybridization  between  races  has  occurred  (Daly  1991;  Del 
Lama  eta/.  1988;  Nunamaker  ef  a/.  1984;  Spivak ef  a/.  1988;  Sylvester  1982, 
1986).  Thus,  supplemental  genetic  characters  subject  to  little  environmental 
modification  have  thus  been  sought.      ^  , 

Mitochondrial  DNA  (mtDNA)  has  been  used  effectively  to  identify  honey 
bee  subspecies  within  populations  in  biogeographic  studies  (Arias,  Soares  & 
Nobrega  1990;  Cornuet  &  Garnery  1991;  Garnery,  Cornuet  &  Solignac  1992; 
Hall  &  Smith  1 991 ;  Meixner,  Sheppard  &  Poklukar  1 993;  Smith  &  Brown  1 988, 
1990;  Smith  1991;  Smith  et  al.  1991).  Honey  bee  mtdna  types  are 
discontinuously  distributed  among  subspecies,  reflecting  the  sustained  temporal 


separation  and  independent  evolution  of  temperate  and  tropical  lineages  (Avise 

et  al.  1987;  Cornuet  &  Garnery  1991).     African  mtDNA  has  been  found  in 

almost  all  feral  neotropical  colonies.  The  small  proportion  of  European  mtDNA 

I 
that  has  persisted  has  remained  largely  confined  to  managed  apiaries  (Hall  and 

Muralidharan  1989;  Hall  &  Smith  1991;  Smith,  Taylor  &  Brown  1989).   These 

findings  supplied  strong  evidence  that  unbroken  matrilines  spreading  as  swarms 

have  been  primarily  responsible  for  the  expansion  of  the  African  bee  population 

.5 

in  the  neotropics  (Hall  &  Muralidharan  1 989;  Hall  &  Smith  1 991 ;  Smith,  Taylor 

I 
&  Brown  1989).   Yet  the  features  of  mtDNA  that  make  these  molecules  good 

genetic   markers   (i.e.    uniparental  inheritance,   no  recombination)   limit  the 

information  that  can  be  obtained  from  their  analysis. 

Genomic  RFLPs  at  several  loci  have  been  identified  that  distinguish  east 

European  from  west  European  and  African  bees,  but  not  west  European  from 

African  bees  (Hall  1 990,  1 992b).  Frequencies  for  the  markers  indicated  limited 

hybridization  of  African  with  European  bees,  but  it  was  not  possible  to  assess     ,,- 

I 
hybridization  with  west  European  bees.   Alleles  at  two  loci,  found  at  different 

frequencies  in  African  and  west  European  bees,  were  found  at  intermediate 

I 
frequencies  in  neotropical  bees,  which  could  have  resulted  from  hybridization 

between  these  two  groups  of  subspecies  (Hall  1992b). 

The  objective  of  the  research  reported  here  was  to  identify  genomic  DNA 

(molecular)  markers  or  characters  that  distinguish  among  the  east  European, 

west  European,  and  African  groups  of  subspecies.  ,    . 
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IDENTIFICATION  AND  GEOGRAPHICAL  DISTRIBUTION 
OF  ALLELES  AT  LOCUS  178 
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Introduction 

Nuclear  DNA  restriction  fragment  length  polymorphisms  (RFLPs)  are 
biparentally  inherited,  codominantly  expressed  markers  shown  to  be  a  valuable 
source  of  genetic  variability  in  honey  bees  (Hall  1986,  1990,  1992b).  Relative 
to  the  biochemical  variability  that  can  be  discerned  in  allozyme  frequencies, 
RFLPs  are  more  abundant  and  do  not  require  gene  expression  for  their 
detection.  East  European-  and  African-specific  RFLP  markers  have  been 
identified,  as  well  as  markers  common  to  African  and  a  west  European  species 
(Hall  1986,  1990,  1992b).  In  neotropical  colonies,  the  east  European  markers 
have  been  found  at  very  low  frequencies  in  areas  where  African  bees  have 
become  established,  indicating  that  there  has  been  limited  paternal  introgression 
from  European  colonies  into  the  African  bee  population  (Hall  1990).  One 
limitation  to  date  in  the  use  of  RFLPs  to  study  ancestry  in  neotropical  bees  has 
been  the  lack  of  additional  markers  specific  for  other  subspecies,  particularly 
A.  m.  mellifera. 

.        Described  in  this  chapter  are  nuclear  DNA  RFLPs,  detected  with  a  single 
genomic  probe,  that  distinguish  bees  of  European  and  African  ancestry.  While 
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it  is  not  possible  to  know  the  genotype  of  each  queen  bee  introduced,  nor 

nnany  details  of  the  events  contributing  to  the  establishment  of  New  World 

honey  bee  populations,  some  reconstruction  can  be  acheived  by  studying 

t 
contemporary  bees  in  the  Old  and  New  World.     Evidence  is  provided  that 

reflects  the  racial  composition  of  Central  and  South  American  bees  prior  to  the 

introduction  of  African  bees,  and  indicates  the  occurrence  of  some,  albeit 

limited,  hybridization  of  European  and  African  bees.  ^  .. , 

I  '  :        '  ;■■■■ 
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Materials  and  Methods 


Sources  of  honev  bees.  Protocols  described  by  Hall  ( 1 986,  1 990)  for  the 
collection,  caste  determination,  and  preservation  of  bees  were  followed  without 
modification.    Drones  were  collected  as  larvae  and  pupae. 

.  South  African  samples,  from  four  locations  in  the  Transvaal,  were 
collected  in  January  1 990  by  HGH.  Brood  samples  from  a  colony  in  So  Frango, 
Brasilia,  were  obtained  in  1 990  by  J.  Maruniak,  University  of  Florida.  Honduran 
samples  included  drones  from  feral  colonies  and  from  managed  colonies 
established  from  feral  swarms,  collected  in  November  1989  by  HGH  and  A. 
Suazo,  Escuela  Agricola  Panamericana,  Honduras.  Samples  from  Tapachula, 
Mexico,  obtained  in  January  1988  by  HGH,  were  from  feral  swarms  captured 
in  bait  hives  maintained  by  the  Mexican  agency  Secretariat  of  Agriculture  and 
Hydrologic  Resources  (SARH),  and  from  two  managed  apiaries.  Sources  of 
USA  drones  included  a  closed  breeding  population  maintained  in  Arizona, 


provided  by  J.  Martin,  G.  Waller,  G.  Loper,  and  E.  Erickson  (Page,  Erickson  & 

I 
Laidlaw  1982;  Severson,  Page  &  Erickson  1986);  managed  colonies  in  Kansas, 

provided  by  0.  Taylor,  University  of  Kansas;  the  University  of  Florida  apiaries; 

and  the  Kona  Queen  Company,  Captain  Hook,  HI.       >  <    ■ 

Electroohoretic  analvsis  of  honev  bee  DNA.  The  cloned  probe  pB  178 
came  from  a  library  of  Ps?l-digested  European  honey  bee  DNA  ligated  into  the 
plasmid  vector  pBR322  (Hall  1986).  This  clone  was  used  as  a  radioactive 
probe  for  detecting  RFLPs  in  honey  bee  DNA  digested  with  restriction 
endonucleases  and  separated  electrophoretically.  Isolation  of  genomic  DNA, 
restriction  endonuclease  digestions,  electrophoresis,  probe  preparation  and 
labeling,  blotting,  and  hybridizations  were  conducted  as  previously  described 
or  cited  (Hall  1986,  1990)  without  further  modification. 

Initial  detection  of  polymorphisms.  In  the  initial  search  for 
polymorphisms,  DNA  was  isolated  from  a  pool  of  ten  workers  from  a  New 
World  European  (USA)  colony  and  from  a  New  World  African  (Costa  Rica) 
colony  (Hall   1986,   1990).     DNA  from  each  pooled  sample  was  digested 


separately  with  nine  restriction  endonucleases,  and  separated  in  2%  agarose 
gels.  The  restriction  fragment  profiles  for  the  European  and  African  samples 
were  compared  following  hybridization  with  cloned  probe  pBI  78.  The  sizes  (in 
kilobase  pairs,  kb)  of  the  restriction  fragments  were  estimated  by  comparison 
to  size  standards  using  a  HI-PAD  digitizer  (Houston  Instruments). 
Polymorphisms    were    detected    in    the    Msp\-    and    Dt/el-treated    samples. 
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Additional  individuals  from  many  locations  in  the  New  and  Old  World  were  then 

I 
analyzed  to  ascertain  the  subspecies  and/or  geographic  distribution  of  the 

polymorphic  fragments. 

The  polymorphisms  generated  by  Msp\  and  Dde\  are  described  here.  For 

the  identification  of  the  Msp\  variants,  402  drone  bees  were  analyzed.  For  the 

identification  of  the  Dde\  variants,  508  drones  were  analyzed,  including  401  of 

those  analyzed  for  the  Msp\  variants.   The  pairwise  combinations  of  the  Msp\ 

and  Dde\  variants  found  in  these  401  drones  are  referred  to  as  alleles. 

Results 

Allelic  nature  of  the  polymorphisms. 

Initially,  polymorphisms  were  detected  with  pB178  following  Msp\  and 
Dde\  digestion  of  DNA  extracted  from  a  sibling  pool  of  New  World  European 
and  African  worker  (female)  honey  bees.  In  the  Msp\  digests,  a  1  .Bkb  fragment 
was  detected  only  in  the  European  sample,  and  a  unique  1.1  kb  fragment  was 
detected  only  in  the  African  sample.  In  the  Dde\  digests,  the  European  sample 
contained  a  unique  1.25kb  fragment,  while  the  African  sample  contained 
unique  3.4kb  and  I.Okb  fragments.  RFLPs  are  codominantly  expressed  in 
diploid  workers,  which  can  result  in  the  comigration  and  masking  of  fragments. 
Therefore,  further  characterization  of  the  Msp\  and  Dde\  RFLPs  was 
accomplished  by  analyzing  drones  (males),  which  are  haploid  offspring, 
parthenogenetically  derived  from  the  queen. 
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The  correspondence  of  the  region  of  the  genome  detected  by  probe 
pBI  78  to  a  contiguous,  heritable  segment  of  a  chromosome,  or  locus  (Fincham 
1983;  Hall  1990,  1992b),  was  demonstrated  by  analyzing  at  least  six  drone 
progeny  from  a  single  colony  (queen).  For  each  pB178/restriction  enzyme 
combination,  no  more  than  two  restriction  fragment  patterns  were  detected, 
which  reflected  the  segregation  of  the  queen's  alleles  (not  shown).  Allelic 
polymorphisms  in  individual  drones  were  manifested  as  different  restriction 
fragment  patterns.  The  different  patterns  generated  by  a  single  restriction 
enzyme  are  referred  to  as  variants.  The  pairwise  combinations  of  the  variants 
produced  by  the  two  enzymes  in  individual  drones  are  referred  to  as  alleles. 

Identification  of  Msd\  variants. 

Most  of  the  36  Msp\  restriction  patterns  or  variants  identified  are  shown 

I 
diagrammatically  in  Figure  1 .   Each  variant  was  composed  of  approximately  8 

to  12  restriction  fragments  between  2.0kb  and  0.3kb.     All  of  the  variants 

included  a  1.4kb  fragment,  one  to  four  fragments  between  1.3kb  and  1.1  kb, 

two  or  three  fragments  of  variable  intensity  between  0.7kb  and  0.5kb,  and  one 

to  four  fragments  smaller  than  0.5kb  (Figures  1,  2,  &  3). 


Msd\  variant  pattern  correlates  with  distribution. 

Variants  with  subsets  of  common  fragments  were  consistently  detected 
in  individuals  from  populations  shown  previously  to  have  similar  genetic 
backgrounds  (Hall  &  Muralidharan  1 989;  Hall  &  Smith  1 991 ).  To  illustrate  this 
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correlation,  variants  consisting  of  patterns  with  common  fragments  and  similar 

•      ■  t 
distributions  have  been  grouped  together.     For  each  population  sampled,  a 

summary  of  the  sample  size,  the  number  of  variants  in  each  Mspl  variant  group, 

and  the  total  number  of  variants  detected  are  given  in  Table  1 . 

-    Variants  Ml 01 -Ml 04  (Figure  2,  first  panel;  see  also  Figure  1)  were 

characterized  by  a  pair  of  fragments  of  approximately  1.2kb  and  0.85kb,  and 

an  intense  signal  at  approximately  0.4kb.    These  variants  differed  from  each 

other  primarily  in  additional  fragments,  between  0.7kb  and  0.5kb.    Variants 

Ml 03  and  Ml 04  were  difficult  to  distinguish  and  will  be  referred  to  together 

as  Ml 03.    Variants  M201  and  M202  (Figure  2,  second  panel)  differed  from 

Ml  01  and  Ml  02  only  in  the  fragments  smaller  than  0.5kb.  Variants  M201  and 

M202  were  distinguished  from  each  other  by  the  fragments  0.7kb  to  0.5kb  in 

size.     The  Ml 00  and  M200  variants  were  found  in  drones  from  the  USA, 

Mexico,  and  Honduras.        "     :       \         )  .,    "      *; 

The  Msp\  variants  M301-M304  (Figure  2,  third  panel;  see  also  Figure  1) 

featured   I.Skb  and   1.2kb  fragments.     Variants  M301,   M303,  and  M304 

differed   very   slightly   from   each   other   in   the   size   of  the   fragments   at 

approximately  1 .4kb  and  1 .2kb  and  in  the  fragments  between  0.7kb  and  0.5kb 

(fragments  not  clear  in  Figure  2  can  be  seen  in  Figure  1).    Variant  M302  was 

distinguished  by  a  0.75kb  fragment  and  by  fragments  smaller  than  0.5kb.    At 

least  one  of  the  M300  variants  was  found  in  each  population. 


-^. 
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Table  1.    Summary  of  populations  examined  for  variants  and  alleles  at  locus  178. 


Sample  origin                             South  Africa 

Honduras 

Mexico 

USA 

Msp\  variant  analysis: 

,t^ 

n:          individuals 

94 

106 

56 

146 

colonies 

23       ^^ 

12 

58 

different  variants: 

19 

11 

8 

Ml 00  total" 

0 

1 

2 

3 

M200  total 

0 

■3  -^  ■■ 

0 

r\  .      3    ^^  ■ 

M300  total 

2 

•r':' 

•   '  ,^  ■ 

2 

M400  total     ■-.■ 

..  3  ■    ■'■•  ... 

./  .3  - . ;. 

''   ■  .5^*^      ■'- 

0 

M500  total 

14 

12 

4 

0 

variants  unique  to 

this  population 

10 

2 

0 

3 

Dde\  variant  analysis: 

"•■.      .    -■-. 

n:         individuals              /■ 

124 

145 

75 

V      164 

colonies 

24 

ZB 

13 

.          59 

different  variants 

19           ' 

.  t?'-^- 

15 

8 

DIOOtotar 

0 

i 

■  :  i:./ 

2 

D200  total 

0 

3 

.- ■■'■!. 

1 

D300  total 

1 

4 

4 

3 

D400  total 

12 

8 

5 

2 

D500  total 

t 

2 

3 

0 

variants  unique  to 

this  population 

10 

3 

1 

"  1 

Allele  identification: 

■  :, 

n:          individuals 

9$ 

106 

56 

146 

colonies 

23 

21 

12 

58 

alleles''  detected 

28 

25 

14 

14 

alleles  unique  to  this 

population 

26 

15 

3 

8 

Total  number  of  different  variants  from  each  variant  group  detected  in  each 
sample  population.         | 

Pairwise  combinations  of  Msp\  and  Dde\  variants  detected  with  pB178  in 
individual  drones. 
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Variants  M401-M405  (Figure  2,  fourth  panel)  were  characterized  by  the 

I 
presence  of  a  1 . 1  kb  fragment  and  three  fragments  smaller  than  0.5kb,  and  by 

^;  .  the  absence  of  the  0.85kb  fragment  common  to  nearly  all  other  variants.  The 

'^:  ...  !■■■       ■  ■  ■  ■   ■,.,:•  ■      ■:■ 

l.lkb  fragment  of  the  M400  group  appeared  to  be  allelic  to  the  0.85kb 
fragment  (Figures  1  &  2).  The  M400  variants  were  distinguished  from  each 
other  on  the  basis  of  restriction  fragments  between  1.4kb  and  l.lkb 
(fragments  smaller  than  0.7kb,  not  clear  in  Figure  2,  can  be  seen  in  Figure  1). 
The  M400  variants  were  found  in  drones  from  South  Africa,  Honduras,  and 


Mexico. 


•^Ci 


^•i^ 


Msp\  variants  M500  as  a  group  (Figures  1  &  3))  were  distinguished  by 

I 
the  presence,  absence,  or  intensity  of  one  or  more  unique  restriction  fragments. 

For  example,  variant  M501  contained  a  unique  1 .5kb  fragment  and  lacked  the 

0.85kb  fragment.   M502  contained  one,  and  M510  and  M51 1  contained  two, 

restriction  fragments  at  approximately  1.3kb;  these  three  differed  from  each 

other  in  the  restriction  fragments  smaller  than  0.7kb.    M503  contained  two 

fragments  near  1.4kb,  similar  to  variant  M403,  but  contained  the  0.85kb 

fragment  that  M403  lacked.    The  M500  variants  were  found  in  drones  from 

South  Africa,  Honduras,  and  Mexico.  >, 

■  ''.■'.  .   ^.  ■.■ 

Identification  of  Dde\  variants.  .  * 

All    33    Dde\    restriction    fragment    variants    identified    are    shown 

diagrammatically  in  Figure  4.  Each  pattern  consisted  of  six  or  seven  restriction 

■  I 
fragments.     Most  variants  exhibited  a  0.75kb  fragment  and  one  to  three 
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fragments  between  1.25kb  and  1.4kb.  The  majority  of  variants  contained  a 
pair  of  fragments  of  variable  size  and  intensity  between  2.0kb  and  2.2kb:  five 
variants  contained  a  2.0kb  fragment,  and  another  fragment  2.3kb  or  larger. 

Many  of  the  Dde\  variants  were  distinguished  by  fragments  that  differed 
slightly  in  intensity  and/or  molecular  weight.  Similar  to  the  Msp\  restriction 
fragment  patterns,  Dde\  variant  groups  were  identified  with  subsets  of  common 
fragments,  and  variants  in  each  group  were  generally  restricted  in  distribution 
to  populations  with  similar  genetic  backgrounds. 


Description  and  distribution  of  groups  of  Dde\  variants. 

Sample  sizes,  numbers  of  variants  in  each  group,  and  total  number  of 
variants  detected  in  each  population  tested  are  given  in  Table  1 . 

Variants  D1 01  and  D1 02  (Figures  4  &  5)  were  characterized  by  a  pair  of 
2. 1  kb  fragments  and  by  1 . 1  kb  and  1  .Okb  fragments.  They  differed  from  each 
other  in  the  sizes  of  fragments  near  2. Okb  and  1.35kb.  These  two  variants 
were  detected  only  in  North  American  drones. 

Variants  D201-D203  (Figures  4  &  5)  were  identified  by  the  presence  of 
a  1.25kb  fragment  and  the  absence  of  the  I.Okb  fragment  detected  in  D101 
and  D102.  Variants  D201-D203  were  distinguished  from  each  other  by  the 
size  and  intensity  of  the  fragments  from  2.1  kb  to  2. Okb  and  1 .35kb  to  1 .25kb 
(Figure  4).  D201  was  found  in  one  USA  colony  and  in  feral  colonies  in  Mexico 
and  Honduras;  D202  and  D203  were  found  in  feral  colonies  in  Honduras. 
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Similar  to  the  D200  variants,  D301-D306  (Figures  4,  5  &  6)  contained 
a  1 .25kb  fragment,  lacked  the  1  .Okb  fragment  detected  In  D1 01  and  D1 02  but 
contained  a  0.75kb  fragment  that  the  D1 00  and  D200  variants  lacked  (Figures 
4,  5  &  6).  D301-D304  were  individually  distinguished  by  the  combinations  of 
fragments  between  2.1kb  and  2. Okb  and  between  1.4kb  and  1.25kb.  D305 
was  further  distinguished  by  a  unique  fragment  of  approximately  1 .8kb.  D306 
contained  a  3.4kb  fragment  and  lacked  a  fragment  near  1.1  kb  found  in  the 
other  D300  variants.  Variants  D301-D305  were  found  in  North  America  and 
Honduras,  while  D306  was  found  only  In  South  Africa. 

The  largest  group,  D400,  consisted  of  fifteen  variants  characterized  by 
I.Okb  and  0.75kb  fragments  (Figures  4  &  7).  The  D400  alleles  differed  from 
each  other  in  the  size  and  intensity  of  fragments  of  approximately  1 .4kb  to 
1.3kb  and  in  the  absence,  or  slight  variation  in  the  size,  of  a  fragment  near 
1.1  kb.  The  majority  of  the  D400  alleles  were  restricted  in  distribution  to  the 
South  African  and  Honduran  populations,  but  two  D400  variants  were  also 
found  in  Mexico.     '  ',  r    ;     {'J     > 

The  D500  variants  were  characterized  by  0.95kb  and  0.75kb  fragments 
(Figures  4  &  5).  Each  variant  in  this  group  differed  in  its  combination  of  two 
fragments  between  2. 1  kb  and  2. Okb,  two  fragments  between  1 .4kb  and  1 .3kb, 
and  a  single  fragment  of  variable  size  near  1  .Ikb.  These  variants  were  found 
in  Mexico,  Honduras,  and  South  Africa. 
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Distribution  of  alleles  at  locus  178  in  Old  and  New  World  drone  honey  bees. 

For  each  population,  the  sample  size,  the  total  number  of  alleles 
detected,  and  the  number  of  unique  alleles  detected  are  given  in  Table  1 .  The 
identification  and  distribution  of  the  65  different  pairwise  combinations  of 
variants,  corresponding  to  alleles,  is  summarized  in  Table  2. 

In  USA  drones,  14  alleles  were  found  that  were  composed  of  Msp\ 

variants  from  the  Ml  00,  M200,  and  M300  groups,  and  Dde\  variants  from  the 

D100,  D200,  D300,  and  D400  groups.     Msp\  variants  Ml 01,  M201,  and 

M203,  and  Dde\  variant  D401 ,  and  a  total  of  eight  alleles  were  unique  to  USA 

drones  (Table  2,  distribution  1). 
•        -  .        I  ■ 

In  the  South  African  drones,  a  total  of  19  Msp\  variants  from  the  M300, 

M400,  and  M500  groups  were  detected;  19  Dde\  variants  from  the  D300, 

D400,  and  D500  groups  were  detected.  Ten  Msp\  and  ten  Dde\  variants  were 

unique  to  the  South  African  samples.     A  total  of  28  different  alleles  were 

identified,  26  of  which  were  not  found  in  any  other  population  (Table  2, 

distribution  8). 

In  samples  from  a  single  colony  in  Brazil,  near  the  site  of  the  introduction 
and  release  of  the  African  bees,  the  Dde\  variants  identified  were  D405  and 
D504,  which  also  had  been  found  in  South  African  drones  (not  shown).  The 
Msp\  variants  in  this  colony,  and  thus  the  alleles,  were  not  determined.      . ; 

In  drones  from  Honduras,  a  total  of  18  Mspl  variants,  17  Dde\  variants, 
and  25  alleles  were  identified.    Alleles  in  these  samples  were  composed  of 
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;    Table  2.   Distribution  of  alleles  detected  at  locus  178  in  Old  and  New  World  drones. 
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Distribution  Key: 


1  =  USA 

2  =  USA  +  Mexico 

3  =  USA  +  Mexico  +  Honduras 

4  =  Mexico 

5  =  Mexico  +  Honduras 

6  =  Honduras 

7  =         ,  Honduras 

8  = 


South  Africa 
South  Africa 
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variants  from  all  five  Msp\  groups  and  from  all  Dde\  variant  groups  except 

D100.  TwoMspl  (M202,  M301)  andtwoDf/el  (D301  and  D201)  variants  were 

common  to  the  Honduranand  USA  samples:  three  Msp I  (I\/I401,  M512,  M514) 

I 
and  five  Dcfe\  variants  (D403-D405,    D407,   D414)   were  common  to  the 

Honduran  and  South  African  drones.    Unique  to  Honduran  drones  were  Msp\ 

variants  M503  and  M518,  and  Dde\  variants  D202,  D203,  and  D305.    Three 

alleles    were    common    to    the    USA    and    Honduras   drones    (M103/D303, 

M301/D201    and  M301/D301;  Table  2,  distribution  3).     Two  alleles  were 

.1 
common    to    Honduran    and    South    African    drones    (M512/D403,     and 

M514/D404;  Table  2,  distribution  7).    Fifteen  of  the  25  alleles  identified  were 

found  only  in  Honduran  drones  (Table  2,  distribution  6).  :  • 

In  bees  from  Mexico,  1 1  Msp\  variants,  1 5  Dde\  variants,  and  a  total  of 

14  alleles  were  identified.    Of  the  Msp\  variant  groups,  only  M200  was  not 

represented,  while  variants  from  all  five  Dde\  variant  groups  were  identified. 

FourMspl  (M102,  M103,  M301,  M302)  and  seven  Dt/el  variants  (D 101,  D102, 

D201,  D301-D303,  D404)  were  common  to  drones  in  Mexico  and  the  USA. 

Three  Msp\  (M405,  M51 1 ,  M51 5)  and  five  Dde\  variants  {D404,  D405,  D408, 

D414,  D501)  were  common  to  drones  in  South  Africa  and  Mexico.    Three 

alleles  were  found  in  both  the  USA  and  Mexico  (M102/D303,  M302/D102,  and 

M301  /D301 ;  Table  2,  distribution  2).  There  were  no  alleles  common  to  South 

African    and    Mexican    drones.       Alleles    M405/D502,     M510/D304,    and 
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M515/D408  (Table  2,  distribution  4),  and  Dde\  variant  D304  were  found  only 

in  this  population.  |  * 

Eight  of  the  65  alleles  contained  exclusive  combinations  of  Msp\  and 

Dde\  variants:    i.e.,  the  Msp\  variant  was  not  associated  with  any  other  Dde\ 

variant.  Fewer  exclusive  combinations  were  found  in  USA/European  bees  than 

in  the  neotropical  and  South  African  bees:   only  M101/D401  was  found  in  the 

USA;  M404/D501    and  M502/D410  were  found  in  Honduras  and  southern 

Mexico;     M513/D504     was     detected     in     Honduras;     and     M304/D503, 

M504/D415,  M505/D507,  and  M519/D402  were  found  in  South  Africa. 

" .    .      •■     \     '^'-      .     V     ■'■.:..•■  1:^  ■'   ■■'  :■■:■■  \;>  -'^:' 
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Common  variants,  uncommon  alleles. 

Several  Msp\  and  Dde\  variants  were  components  of  numerous  alleles 
and  were  found  in  more  than  one  population  (Table  2).  The  combination  of 
variants  comprising  each  allele  was  indicative  of  the  ancestry  of  the  sample. 
For  example,  Dde\  variant  D405  was  found  in  seven  alleles  in  South  Africa  and 
the  neotropics.  The  majority  of  samples  with  D405  were  from  South  Africa, 
and  the  associated  Msp\  variants  in  these  individuals  (M400s  and  M500s)  were 
only  found  in  Old  and  New  World  colonies  previously  shown  to  have  African 
mtDNA  (Hall  &  Muralidharan  1989;  Hall  &  Smith  1991).   In  North  America  and 
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Honduras,  Dde\  variant  D303,  found  in  three  alleles,  was  associated  with  Msp\ 
variants  M102,  M103,  and  M201,  indicating  European  origin. 

The  presence  of  Dde\  variant  D404  in  the  USA,  Honduras,  and  South 
Africa  was  an  exception  to  the  correlation  of  fragment  pattern  with  geographic 
origin.  The  alleles  with  D404  were  differentiated  by  Msp\  variants,  which  were 
consistent  with  the  sample  origin  and  what  was  known  of  its  ancestry:  M302 
(USA),  M51 4  (Honduras  and  South  Africa),  and  M51 5  (South  Africa)  (Table  2). 

Variants  in  the  M300  group  were  found  in  all  populations  examined, 

although  individually  they  appeared  to  be  limited  in  distribution.  Variant  M302, 

which  contained  a  0.75kb  fragment  not  present  in  M301,  M303,  and  M304, 

I    ■■        ■ 
was  found  only  in  the  USA  and  Mexico,  with  Dde\  variants  D101,  D102,  and 

D404.    Variant  M301  was  a  component  of  two  different  alleles  in  the  USA, 

with  Dde\  variants  D201  and  D301 .  Variants  M303  and  M304,  very  similar  to 

M301  but  found  only  in  South  Africa,  were  associated  with  Dde\  variants  D405 

and  D503  (Figures  1  and  2).  •  f  *  -»•   ^. 

In  neotropical  samples,  because  of  the  small  difference  in  fragment  sizes 

that  distinguished  M301,  M303,  and  M304,  the  African  or  European  origin  of 

i 

the  alleles  containing  the  Mspl-1.8kb  fragment  was  confirmed  by  the  Dde\ 
variant  with  which  they  were  associated.  As  in  the  USA,  M301  was  found 
with  D201  and  D301  in  drones  in  Mexico  and  Honduras,  indicating  that  the 
samples  were  European.  In  addition,  variants  D202,  D203,  and  D305,  found 
only  in  drones  from  Honduras,  had  fragments  common  to  variants  found  in  USA 
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drones.      Given  that  the  fragment  patterns  of  these  Dde\   variants   were 

consistent  with  other  patterns  detected  in  the  USA,  these  three  variants  were 

assumed  to  indicate  European  ancestry.  Neotropical  drones  in  which  the  Msp\- 

1.8l<b  fragment  was  detected  were  concluded  to  be  of  European  ancestry  on 

the  basis  of  the  Dde\  variants  (Table  2).  - 
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Greater  diversitv  at  locus  1 78  in  tropical  vs.  temperate  bees. 

USA  bees  demonstrated  the  least  variability,  having  the  smallest  number 

of  variants  and  alleles  found,  of  the  four  populations  analyzed.    The  greatest 

diversity  at  locus  1 78  was  found  in  the  South  African  drones:  twice  as  many 

Mspl  and  Dde\  variants  comprising  twice  as  many  alleles  were  found  in  half  as 

many  colonies,  compared  to  the  USA. 

The  neotropical  samples  were  also  more  varied  than  the  USA  samples 

I  ■ 

and  were  collected  from  fewer  colonies.    The  Honduran  samples  came  from 

colonies  established  from  feral  swarms  several  years  after  African  bees  had 

arrived.   The  detection  of  both  European  and  African  variants  and  alleles  was 

expected.    However,  there  was  little  indication  of  European  ancestry  in  the 

Honduran  bees;  the  majority  of  alleles  common  to  the  USA  and  Honduras 

contained  the  Mspl   variant  M301,   and  only  two   other  alleles  contained 

USA/European-type  variants  (Table  2,  distribution  3).      More  variants  were 

common  to  Honduran  and  South  African  drones.     Like  the  South  African 

drones,  the  majority  of  the  Mspl  variants  in  the  Honduran  drones  belonged  to 

the  M400  and  M500  group,  and  most  of  the  Dde\  variants  detected  were  in  the 
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D400  group.  This  finding  was  consistent  with  the  feral  origin  of  the  colonies, 

in  which  African  mtDNA  was  found  (Hall  &  Smith  1991). 

Some  of  the  drones  from  southern  Mexico  were  collected  from  managed 

colonies  established  from  feral  swarms,  while  the  majority  were  collected  from 

feral    swarms    (Hall    &    Muralidharan    1989).      It   was   expected   that   both 

USA/European  and  South  African  markers  would  be  found.    The  number  of 

variants  common  to  Mexico  and  South  Africa  was  far  greater  than  those 

common  to  the  USA  and  Mexico,  although  there  were  no  alleles  common  to  the 

I 
Mexican  and  South  African  sample  populations.  The  variants  and  alleles  unique 

to  the  Mexican  samples  resembled  those  in  the  South  African  samples  and, 

apparently,  represented  the  recent  establishment  of  the  feral  African  population. 

As  in  Honduras,  most  of  the  European-type  variants  and  alleles  common  to  the 

USA  and  Mexican  populations  were  in  the  M300  group  (Table  2,  distributions 

2  and  3). 

The  greater  number  of  variants  and  alleles  identified  in  the  neotropical 
samples  relative  to  the  USA  reflected  the  presence  of  African  as  well  as 
European  variants.  Variants  unique  to  the  neotropical  drones,  but  similar  to 
those  found  in  USA  bees,  probably  reflect  regional  differences  in  European 
races.  Five  alleles  detected  only  in  Honduras  and  Mexico  (Table  2,  distribution 
5)  were  composed  of  variants  that  either  were  found,  or  were  of  the  same 
variant-type  as  those  found,  in  the  South  African  drones  (Table  2,  distribution 
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8).   The  detection  of  these  alleles  suggests  that  greater  variation  exists  in  the 
parent  South  African  population  than  has  been  revealed  in  this  study. 

Alleles  were  identified  in  Mexico  (M510/D304)  and  Honduras 
(M202/D414,  and  M510/D302)  that  were,  according  to  the  assignment  of 
variants  to  groups  on  the  basis  of  fragment  pattern  and  distribution,  composed 
of  both  European  and  African  variants.  The  distributions  determined  in  this 
study  may  not  be  absolute,  given  the  number  of  variants  and  alleles  found,  the 

size  of  each  sample,  and  the  world-wide,  human-assisted  redistribution  of  bees. 

I 
Investigation  with  additional  samples  will  likely  be  needed  for  a  thorough 

classification  of  variants  at  this  highly  polymorphic  locus.    However,  alleles 

composed  of  African  and  European  variants  may  be  the  result  of  recombination 

within  the  178  locus. 


,  ^ 


Greater  allele  diversity  has  been  detected  in  South  African  compared  to 
USA  honey  bees  at  another  genomic  locus  (Hall  1992b).  At  least  14  variants 
were  detected  in  South  African  bees,  two  of  which  were  detected  in  European 
and  USA  bees.  In  Honduras  nine  variants  were  detected:  the  two  common  to 
European  and  USA  samples,  as  well  as  seven  characteristic  of  South  African 
bees.  A  larger  number  of  mtDNA  size  classes  have  been  found  in  African  bees 
than  in  bees  of  European  ancestry  (Hall  &  Smith  1991;  Smith  et.  al.  1991). 
Greater  diversity  has  been  revealed  in  tropical  relative  to  temperate  populations 
of  Drosophila  melanogaster  (David  &  Capy  1988;  Hale  &  Singh  1987), 
Drosophila  simulans  (Hyytia  et  al.    1985),  Ceratitis  capitata  (Gasperi  et  al. 
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1991),  and  Limulus  polyphemus  (Saunders,  Kessler  &  Avise  1986).     This 

diversity  has  been  attributed  to  the  capacity  of  tropical  populations  to  maintain 

large  population  sizes  over  a  long  period  of  time  (David  &  Capy  1 988)  and  may 

be  consistent  with  the  large  reproductive  capacity  and  population  sizes  of 

tropical  bees  relative  to  temperate  bees  (Winston,  Taylor  &  Otis  1983). 
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ALLELE  FREQUENCIES  AT  LOCUS  178  REVEAL  HYBRIDIZATION 
OF  EUROPEAN  AND  AFRICAN  BEES  IN  THE  NEOTROPICS 


Introduction 
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Frequencies  are  reported  here  for  the  individual  variants,  for  the  groups 
of  variants  with  connmon  fragments  and  distributions,  and  for  the  alleles  found 
at  locus  1 78  in  Old  and  New  World  populations  of  honey  bees.  Unlike  previous 
nuclear  DNA  markers  (Hall  1 990,  1 992b),  variants  and  alleles  detected  at  locus 
1 78  appeared  to  be  specific  to  east  European,  west  European,  or  African  bees. 
Ancestry  in  New  World  bees  was  inferred  from  the  frequencies  for  the  variants, 
variant  groups,  and  alleles.  The  results  reveal  a  level  and  specificity  of  African- 
European  hybridization  not  observable  in  previous  DNA  studies. 


Materials  and  Methods 
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Procedures  for  the  isolation  and  electrophoretic  analysis  of  honey  bee 
genomic  DNA,  identification  of  the  source  and  size  of  probe  pBI  78  used  for  the 
detection  of  RFLPs,  and  the  details  of  the  initial  detection  of  polymorphisms, 
were  given  in  the  preceeding  chapter.  RFLPs  were  initially  detected  in  pooled 
sibling  samples  of  worker  bees.   For  determining  the  frequencies  for  the  Msp\ 
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and  Dde\  variants  detected  with  pB178,  individual  Old  and  New  World 
European  and  African  bees  were  analyzed. 

Sources  of  honey  bees.  Frequencies  of  variants  and  alleles  in  drones 
were  determined  for  the  samples  used  for  the  variant  and  allele  identification. 
Worker  larvae  and  pupae  were  collected  from  these  same  colonies  in  South 
Africa,  Honduras,  and  southern  Mexico. 

Adults  workers  from  Venezuela,  collected  between  1986  and  1988, 
were  provided  by  R.  Hellmich  Jr.,  J.  Villa,  A.  Collins,  and  T.  Rinderer,  USDA- 
ARS,  Baton  Rouge.  Costa  Rican  samples  were  obtained  in  May  1989  from 
apiaries  at  Cerro  de  la  Muerte  and  San  Isidro  del  General,  by  HGH  with  the  help 
of  H.  Arce  and  R.  Dormond,  National  University,  Heredia,  Costa  Rica.  Five 
colonies,  maintained  at  2200m  to  test  resistance  to  African  introgression  at  a 
higher  elevation  had  European  mtDNA.  Three  colonies  maintained  at  700m  had 
African  mtDNA.  Samples  from  southern  Mexico,  near  Tapachula,  were 
obtained  in  January  1988  by  HGH  from  two  managed  apiaries  and  from  feral 
swarms  captured  in  bait  hives  maintained  by  the  Mexican  agency  Secretariat 
of  Agriculture  and  Hydrologic  Resources  (SARH).  Feral  worker  samples  from 
northern  Mexico,  collected  prior  to  the  arrival  of  African  bees,  were  provided 
by  W.  Rubink  and  A.  Collins,  USDA-ARS,  Weslaco.  Some  of  the  USA  drones, 
and  all  of  the  USA  workers,  were  from  a  closed  breeding  colony  in  Arizona, 
provided  by  J.  Martin,  G.  Waller,  G.  Loper,  and  E.  Erickson  (Page,  Erickson  & 
Laidlaw  1982;  Severson,  Page  &  Laidlaw  1986).    Sources  for  additional  USA 


38 
drones  were  reported  in  the  previous  cfiapter.  Workers  fronn  Europe  (larvae  and 

pupae)  were  provided  by  B.  Vaissiere,  Texas  A&IVI,  and  J.-M.  Cornuet,  INRA 

Monfavet,  France.   These  included  samples  of  west  European  A.  m.  mellifera 

(8  colonies),  and  east  European  A.   m.   ligustica  (3  colonies)  and     A.   m. 

caucasica  (2  colonies). 

For  determining  the  frequencies  of  the  Msp\  variants,  a  total  of  862  bees 

were  analyzed  (402  drones,  460  workers),  representing  128  Old  and  New 

World  colonies.   For  determining  the  Dde\  variant  frequencies,  the  402  drones 

analyzed  for  the  Msp\  variant  frequencies  and  an  additional  106  drones  were 

examined,  from  a  total  of  1 22  colonies.  Allele  frequencies  were  determined  for 

I 
401  drones  from  1 14  colonies,  for  which  both  the  Msp\  and  Dde\  variants  had 

been  identified. 

Frequencies  of  variants  and  alleles.    In  haploid  drones,  frequencies  for 

Msp\  and  Dde\  variants  and  alleles  were  calculated  as  a  fraction  of  the  total 

number  of  variants  or  alleles  detected.    The  numbers  of  variants  and  alleles 

detected  were  determined  as  described  by  Hall  (1992b),  using  the  data  for 

drones.  At  least  two  drones  were  analyzed  per  colony.  When  fewer  than  six 

■  I  ^■ 

drones  were  analyzed  per  colony  and  only  one  variant  or  allele  was  found,  the 

variant  or  allele  was  counted  once.  When  more  than  six  drones  were  analyzed 

and  one  variant  or  allele  was  detected,  it  was  counted  twice  (greater  than  95% 

probability  that  the  queen  was  homozygous).    When  two    variants  or  alleles 


were  found  (heterozygous  queen),  each  was  counted  once.    Frequencies  for 

groups  of  variants  represent  the  sum  of  the  individual  variants  within  a  group. 

Drones  from  the  same  colony  represent  only  the  variants  and  alleles  of 

I       ■ 
a  single  individual,  the  queen.    Frequencies  for  variants  identified  from  the 

analysis  of  worker  restriction  fragment  patterns  reflect  more  of  the  variation 
within  a  local  population,  due  to  the  presence  of  multiple  patrilines  within 
colonies.  i 

In  diploid  worker  bees  (females),  RFLPs  are  codominantly  expressed. 
When  two  restriction  fragment  variants  contain  common  fragments,  fragment 
superimposition  results  as  a  consequence  of  comigration.  Thus  for  workers, 
frequencies  were  determined  for  groups  of  Msp\  variants,  since  more  than  one 
variant  from  the  same  group  could  account  for  the  fragment  patterns  seen 
(Figure  8),  and  since  the  variant  group  was  correlated  to  the  ancestry  of  the 
sample  (see  previous  chapter).  If  the  identity  of  the  variant  group  to  which  a 
particular  fragment  belonged  was  uncertain,  the  frequency  for  the  group  was 
reported  as  a  range  (Table  4).  '       .  .     •'  ''        ';r 

In  Dde\  digests  of  the  majority  of  worker  samples,  fragment  comigration 

t 
precluded  the  identification  of  variants  (Figure  1 1 )  and  the  determination  of  the 

frequencies  of  both  individual  variants  and  of  variant  groups.    Consequently, 

frequencies  of  Dcfe\  variants,  and  therefore  of  alleles,  were  determined  only  for 

those  populations  for  which  drones  were  available  (South  Africa,  Honduras, 

southern  Mexico,  and  the  USA;  Tables  5  &  6). 
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Results 


""■"  Frequencies  determined  in  drones  for  the  individual  Mspl  restriction 
fragment  patterns  (variants),  and  totals  for  the  five  groups  of  Msp\  variants 
with  similar  fragment  patterns  and  distributions  are  presented  in  Table  3.  The 
fragments  characteristic  of  each  group  of  Msp\  variants  are  shown  in  Figure  8. 

Frequencies  for  groups  of  Msp\  variants  detected  in  workers  are 
summarized  in  Table  4.  Examples  of  Msp\  restriction  fragment  patterns  are 
shown  in  Figure  9.  The  pair  of  Msp\  variant  groups  detected  in  a  worker  will 
be  referred  to  as  the  Mspl  genotype.  Three  worker  restriction  fragment 
patterns  shown  in  Figure  10  demonstrate  how  the  comigration  of  fragments 
confounded  the  identification  of  individual  Msp\  variants  in  diploids. 

Frequencies  for  Dde\  variants  detected  in  drones  are  summarized  in  Table 
5,  and  the  fragments  characteristic  of  each  of  the  five  Dde\  variant  groups  are 
shown  in  Figure  1 1 .  Dde\  fragment  patterns  were  examined  in  selected 
workers  for  evidence  of  fragments  characteristic  of  particular  variant  groups 
(not  shown).  '      '  ,        :      v 

The  distribution  and  frequencies  for  the  alleles,  the  pairwise  combinations 
of  the  Mspl  and  Dde\  variants,  detected  in  drones  are  summarized  in  Table  6. 
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European  bees. 

A  limited  number  of  workers  were  available  from  European  colonies  of 
known  ancestry.    Thus,  the  variants  represent  a  subset  of  the  variability  in 
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Table  3. 

Frequency  [%]  of  detection  of  Msp\  variants  at  locus  178 

in 

drones. 

Variant 

South  Africa         Honduras              Mexico 

USA 

M101 
M102 
M103 

Z                       19 

3 
13 
51 

M201 

11 

M202 

.   ..        .,;-         1 

:      ^ 

3 

M203 

'     -                      1 

1 

M301 

1 

26 

24 

9 

M302 

-.-.  _        1 

10             ., 

9 

M303 

3 

M304 

3 

M401 

3 

3 

M402 

5 

K''  \.   ' 

M403 

1 

..V-       t 

10 

M404 

8 

5 

M405 

5 

5       '^' 

M501 

3 

M502 

5 

•   ■  ■■■  -S     ■.^--   > 

M503 

3 

M504 

3 

.'"■■"          ■',-'■■■',■- 

M505 

3 

M506 

-        V         5 

.       ■        .           A         ('■■■' 

M507 

■        i          5 

M508 

27 

3 

M509               •  V.  y  ' 

3 

,815 

M510                    ' 

:^    ■              / 

'  •   v  6 

•        •■■       '■;     „5            , 

M511 

3 

8 

5    ■■. 

*    ' 

M512               r 

.rV         ""l 

10 

,^,i\ri 

M513 

3 

•^  ■  X'''!.    -,V  .^ 

M514 

5 

5 

■  ./"•■-■i  v-f, 

'!■':  •<>.    r 

M515 

3 

3 

-  ■  -5    ■■  ^.^  ■ 

•■  ,v         -         -    ■ 

M516"' 

5 

3 

;■-!-  4' 

M517"' 

3 

""':>';.      - 

M5  IS- 

3 

-/^^;  .• 

MS  19- 

3 

Ml 00  total" 

0 

3 

24 

67 

M200  total 

0 

3 

0 

14 

M300  total 

5 

26 

33         ;, 

18 

M400  total 

14 

18 

19        "vA 

0 

M500  total 

81 

51 

24 

0 

n:    individuals  94 

colonies  23 

variants  detected''       37/44 


106 
21 

39/42 


56 

12 

21/24 


146 

58 

76/112 


^^;^- 


ns 
a 
b 


Not  shown  in  a  figure. 

Sum  of  the  frequencies  of  all  variants  detected  in  the  respective  variant  groups. 
Ratio  of  the  number  of  variants  counted,  to  the  total  number  of  variants 
possible  =  I[no.  drones/colony]. 
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Figure  8.  Restriction  fragment  patterns  of  Msp\  variants 
detected  at  locus  1 78.  Restriction  fragments  detected  in 
all  variants  within  a  group  and  characteristic  of  each 
group  are  indicated  by  (_),  and  variable  fragments 
detected  in  at  least  one  but  not  all  variants  within  each 
group  of  Msp\  variants  are  designated  by  (D). 
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Size 

(kb) 

M403 
M503 

#2 

M301 
M202 

M301 
M509 

M304 
M512 

#3 

M404 
M201 

M404 
M508 

M401 
M502 

2.0kb 
1.8kb 

1.4kb 

1.2kb 
l.lkb 

I.Okb 

O.BBkb 
0.75kb 

0.6kb 
0.4kb 

i 

• 

III          III                      1  1     1 1 
II       II        III        1                III 

m 

II  1            III              1    1     1     1 
1             II                       1              III              1    1     1     1 

III  1              III              1    1     1     1 
1              II                       1            III              1    1     1     1 

1         II                       1                  II              1    1     1     1 
1             II                      1           III             1    1     1     1 

1 

III      III             1 1 II   1 
1  1  1  1      III     1          II    1 

III      III             1  III   1 
1  1  1  1      III     1          II    1 

III      III             III 
1  1  1  1      III      1           II 

Figure  10.  Msp\  patterns  from  three  workers,  demonstrating  possible  variant  pairs. 
#1 :  an  example  of  an  unequivocal  identification  of  the  variants  which  composed  the 
pattern  of  fragments  detected.  #2  and  #3:  examples  demonstrating  the  difficulty 
encountered  in  identifying  variants  in  workers  from  neotropical  areas  where  African 
bees  were  established.  Fragments  in  these  patterns  could  be  assigned  to  variants  in 
more  than  one  variant  group. 
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Table  5.    Frequency  [%]  of  detection  of  Dde\  variants  at  locus  178  in  drones. 


iy» 


Variant 

South  Africa 

Honduras 

So.  Mexico 

USA 

D101 

4 

1 

D102 

■  • 

4 

6 

D201 

2 

4 

7 

D202 

1 

2 

D203 

2 

D301 

13 

17 

1 

D302 

:     ■.•;'' 

4 

4 

28 

D303 

2 

25 

53 

D304 

4 

D305 

"  Ji:-':! 

,  2 , :, 

,         >- 

D306 

' 

D401 

2 

D402 

2 

■       -"'.    4   -    •'.' 

D403 

13 

19 

D404          ... 

2 

4 

1 

D405 

22 

6-8 

g 

D406 

2 

4 

D407 

11 

8-11 

>v^. 

D408 

9 

•4     ...-. 

D409 

2 

.    .^.^ 

D410 

4 

4 

D411 

2 

D412 

2 

D413 

2 

•  .■■.,. 

D414                        ^ 

11 

13 

4 

D415 

2 

D501 

2 

6 

4 

D502 

/    ,.■         2 

"■"'.'"■ 

^H:'- 

D503 

•...  .V.         2 

D504 

t>. 

2 

4 

D505 

4-7 

■  .,'' 

D506 

0-2 

D507 

2 

D1 00  total" 

0 

0 

8 

7 

D200  total 

0 

6 

4 

7 

D300  total 

2 

21 

44 

82 

D400  total 

^.r        82 

64 

25 

4 

D500  total 

16 

8 

12 

0 

n:  individuals 

124 

145 

75 

164 

colonies 

24 

26 

13 

59 

variants  detected'' 

45/48 

47/52 

24/26 

83/118 

«  Total  number  of  different  variants  from  each  group  detected  in  each  population. 

b  Number  of  variants  counted/total  number  of  variants  possible  (two  per  queen 

or  colony). 
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Figure  1 1 .  Restriction  fragment  patterns  of  groups  of 
Dde\  variants  detected  at  the  locus  corresponding  to 
probe  pB178  in  the  honey  bee  genome.  Restriction 
fragments  detected  in  all  variants  within  a  group  and 
characteristic  of  each  group  are  indicated  by  (-),  and 
variable  fragments  detected  in  at  least  one  but  not  all 
variants  within  each  group  of  Dde\  variants  are  designated 
by  (D).  ,     ,.  . 
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these  east  and  west  European  populations.    There  was  little  variation  in  the 

■.-  -     .  ,-        '-      -f;'  ' 

restriction  fragment  patterns  within  each  European  sample  population. 

Examples  of  Msp\  fragment  patterns  detected  in  A.  m.  ligustica  workers 
are  shown  in  Figure  9,  lanes  1-3.  The  fragments  detected  corresponded  to 
Ml  00  and  M200  variants,  with  estimated  group  frequencies  of  81  %  and  19%, 
respectively  (Table  4).  The  majority  of  A.  m.  ligustica  workers  appeared  to 
have  two  Ml 00  variants  (Figure  9,  lane  1),  while  the  A.  m.  caucasica  workers 
examined  all  appeared  to  have  two  M200  variants  (Figure  9,  lane  3).  The  Dcle\ 
fragment  patterns  for  the  A.  m.  ligustica  and  A.  m.  caucasica  workers 
examined  (not  shown)  contained  1.25kb  and  0.75kb  fragments  characteristic 
of  variants  in  the  D300  group  (Figure  11).  It  is  possible  that  D200  variants 
were  present  but  could  not  be  detected  due  to  comigration  of  D300  fragments 
(see  Figures  4  &  11).  '^'s^'-k 

Examples  of  Msp\  fragment  patterns  detected  in  A.  m.  mellifera  workers 
are  shown  in  Figure  9,  lanes  4-7.    All  the  A.  m.  mellifera  workers  examined 
•^U;  appeared  to  have  at  least  one  M300  variant.    The  frequency  for  the  M300 

variant  group  was  94%  (Table  4).  The  majority  of  the  workers  appeared  to 
have  two  M300  variants  (Figure  9,  lanes  6  &  7).  The  remaining  A.  m.  mellifera 
, :,/.  workers  appeared  to  have  an  M300  variant  together  with  an  Ml 00  or  M200 
variant  (Figure  9,  lanes  4  &  5,  respectively).  All  Ddel  fragment  patterns 
detected  in  the  A.  m.  mellifera  workers  contained  a  1.25kb  fragment 
characteristic  of  the  D200  and  D300  groups  (not  shown).     The  majority 


(-60%)  of  the  patterns  in  the  A.  m.  mellifera  samples  contained  1.25kb  and 
0.75kb  fragments  (the  1.25l<b  fragment  was  not  detected  in  South  African 
drones),  and  lacked  a  fragment  at  approximately  1  .Okb  detected  in  the  majority 
I8.  of  South  African  drones  (Figures  4,  6,  &  7).  Ihe  A.  m.  mellifera  Dde\  ^atietns 

,  lacking  the  0.75kb  fragment  ( ~  25%)  appeared  to  be  composed  of  one  or  two 

D200  variants.  In  a  few  (-5%)  workers  the  1.25kb  and  1  .Okb  fragments  were 
present,  but  the  0.75kb  fragment  was  absent,  indicative  of  D100  and  D200 

»■.: 

':  variants,  which  were  not  detected  in  South  African  drones.  A  combination  of 

D100,  D200,  and  D300  variants  was  indicated  in  the  workers  (10%)  in  which 

■  the  1.25kb,  I.Okb,  and  0.75kb  fragments  were  detected.    No  evidence  was 

found  in  the  A.  m.  mellifera  workers  for  Dde\  patterns  detected  in  South  African 

<  .'         drones  (Figures  4  &  1 1). 

African  bees.  '  . 

;'  In   South   African  drones,  the  Msp\   variant  group  with  the   highest 

frequency  (81%)  was  the  M500  group  (Table  3).  Variants  M508  and  M512 
had  the  highest  individual  frequencies  (25%  and  1 1%,  respectively).  Eleven 
other  M500  variants  were  found  at  lower  frequencies,  each  in  one  or  two 
drones.  The  M400  group  frequency  was  14%:  variants  M401,  M402,  and 
M405,  were  each  detected  at  a  low  frequency.  The  frequency  for  the  M300 
group  was  6%:  variants  M303  and  M304  were  each  detected  in  a  single  drone 
from  colonies  in  different  locations  in  the  Transvaal  (Table  3). 
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53 
Examples  of  Msp\  fragment  patterns  detected  in  South  African  workers 

are  shown  In  Figure  9,  lanes  12,   14,   16,  and  17.     The  estimated  group 

frequency  for  the  M500  variants  In  South  African  workers  was  92%  (Table  4): 

the  majority  of  these  workers  appeared  to  carry  only  M500  variants  (Figure  9, 

lane   17).      The   M400   and   M300   group  frequencies  were   5%   and   3%, 

respectively.   Since  no  M200  variants  were  detected  In  South  African  drones, 

it  was  assumed  that  M500  variants  predominated  In  workers  as  well  (Table  4: 

M200  variants  would  be  difficult  to  detect  In  A.  m.  scutellata  workers  due  to 

the  comlgratlon  of  fragments).  ;. 

In  South  African  drones,  the  majority  of  the  Dde\  variants  detected  were 
In  the  D400  group  (Table  5).  Variants  D405  and  D403  were  found  at  the 
highest  individual  frequencies  (22%  and  13%,  respectively).  Only  one  D300 
variant  was  detected  (D306).  The  remaining  variants  detected  were  from  the 
D500  group.  The  ambiguity  In  the  frequencies  for  D505  and  D506  was  due  to 
the  similarity  of  these  patterns,  which  were  difficult  to  distinguish  If  the 
samples  were  not  adjacent  on  a  blot. 

Frequencies  were  low  for  the  many  alleles  found  In  South  African  drones 
(Table  6).  Alleles  l\/I508/D408  and  M512/D403  were  detected  at  the  highest 
frequencies  (1 1  %  and  8%,  respectively):  two-thirds  of  the  alleles  were  each 
detected  In  a  single  drone  (equivalent  to  one  allele  In  the  queen  from  a  single 
colony,  or  one  In  forty-six).  _ 
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North  American  bees.  *  '       :■      . 

United  States.  In  USA  drones,  the  variant  group  with  the  highest 
frequency  was  IVI100  (67%),  and  the  Msp\  variant  detected  at  the  highest 
frequency  was  M103  (51%;  Table  3).  Variant  M201  was  detected  at  the 
highest  frequency  (1 1  %)  in  the  M200  group.  Variants  M301  and  M302  were 
each  detected  at  a  frequency  of  9%. 

In  USA  workers,  Msp\  fragment  patterns  characteristic  of  the  Ml 00, 
M200,  and  M300  variants  were  detected  at  group  frequencies  of  80%,  5%, 
and  14%,  respectively  (Table  4).  The  majority  of  workers  appeared  to  have 
two  Ml 00  variants  (Figure  9,  lane  1).  All  other  Msp\  fragment  patterns  in  USA 
workers  corresponded  to  heterozygous  genotypes  (Figure  9,  lanes  2,  4,  and  5, 
respectively). 

In  drones  from  the  USA,  Dde\  variants  D302  and  D303,  and  the  D300 
group,  were  detected  at  the  highest  individual  and  group  frequencies, 
respectively  (Table  5).  Frequencies  were  low  for  all  the  other  variants  detected. 

Allele  M103/D303  was  detected  at  the  highest  frequency  (39%)  in 
drones  from  the  USA,  the  highest  frequency  of  all  alleles  in  all  Old  and  New 
World  samples  examined.  Frequencies  for  the  other  alleles  were  relatively  low 
(Table  6).  Half  as  many  alleles,  in  twice  the  number  of  colonies,  were  detected 
in  USA  drones  compared  to  South  African  drones. 

Northern  Mexico.  Msp\  fragment  patterns  detected  in  feral  workers  from 
northern  Mexico  are  shown  in  Figure  9,  in  lanes  2,  5,  6,  and  7.    The  variants 
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detected  were  from  the  Ml  00,  M200,  and  M300  groups,  at  group  frequencies 

of  9%,  27%,  and  64%,  respectively.    Only  in  the  A.  m.  mellifera  workers  from 

France  was  the  M300  group  frequency  higher.         ^    -r  "'    . 

'     -  .  -i".-  ■ 
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Neotropical  bees.  ty    '  J^;  . 

Southern  Mexico.  Variants  from  four  of  the  five  Msp\  variant  groups 
were  detected  in  drones  from  southern  Mexico,  collected  after  African  bees 
were  established.  The  group  with  the  highest  frequency  was  M300  (33%). 
The  variant  with  the  highest  individual  frequency  was  M301  (24%).  As  in  USA 
drones,  M302  was  the  only  other  M300  variant  detected.  Frequencies  for  the 
Ml  00,  M400,  and  M500  variant  groups  were  comparable  to  each  other  (Table 
3).  Of  the  Ml 00  variants.  Ml 03  was  detected  at  the  highest  frequency 
(19%),  as  in  the  USA  drones.  Frequencies  were  low  for  the  individual  M400 
;^-3^;  and  M500  variants  detected  (Table  3).      ..      ;  «  r 

Examples  of  Msp\  fragment  patterns  detected  in  workers  from  southern 
Mexico  are  shown  in  Figure  9,  in  lane  2  and  in  lanes  4-17.  In  these  samples, 
fragments  characteristic  of  variants  from  all  five  Msp\  variant  groups  were 
found  (Table  4).  Frequencies  for  the  M300,  M400,  and  M500  variant  groups 
were  32%,  21  %,  and  30%,  respectively.  The  Ml 00  and  M200  variant  groups 
were  detected  at  the  lowest  frequencies. 

Variants  from  all  five  Dde\  groups  were  detected  in  drones  from  southern 
Mexico.  The  Dde\  variants  detected  at  the  highest  frequencies  were  D303  and 
D301 .    Frequencies  were  low  for  all  other  variants  detected  (Table  5). 
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Allele  frequencies  in  drones  from  southern  Mexico  were,  in  general, 

intermediate  to  the  frequencies  in  the  USA  and  South  African  drones.    The 

notable  exception  was  allele  M301/D301  (19%),  which  was  absent  in  South 

African  drones  and  detected  in  a  single  colony  from  the  USA  sample  collection 

(a  colony  from  Florida).  ^ 

Honduras.  Variants  from  each  Msp\  group  were  detected  in  drones  from 
Honduras  (Table  3).  The  majority  of  variants  detected  were  in  the  M500  group 
(Table  3).  With  the  exception  of  variant  M51 2,  detected  at  a  frequency  of  1 1- 
14%,  the  individual  M500  variants  were  found  at  low  frequencies.  [An 
ambiguity  in  the  frequencies  of  variants  M512  and  M516  resulted  from  the 
similarity  in  the  sizes  and  intensities  of  fragments  in  these  two  variants:  unless 
the  samples  bearing  these  variants  were  adjacent  in  Southern  blots,  they  were 
difficult  to  distinguish.]  Frequencies  of  variants  in  the  M400  group  (M401 ,  3%; 
M403,  8%;  M404,  8%)  were  comparable  to  the  frequencies  at  which  these 
variants  had  been  detected  in  South  Africa  and  southern  Mexico  (Table  3).  As 
in  southern  Mexico  and  the  USA,  only  M301  and  M302  were  detected  from  the 
M300  group.   Variants  Ml 03  and  M202  were  detected  in  single  drones. 

It  appeared  that  all  five  Msp\  variant  groups  were  represented  in  workers 
from  Honduras  (Figure  9,  lanes  4,  6,  8,  11,  12,  14-17;  see  Table  4).  As  in 
drones,  the  variant  group  detected  at  the  highest  frequency  was  M500  (39%). 
The  frequency  for  the  M400  group,  27%,  was  the  highest  of  all  the  New  World 
populations  tested.    The  M300  group  frequencies  (26%  in  drones,  23%  in 
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workers)  were  comparable  to  other  neotropical  colonies  examined.  The  Ml  00 

and  M200  groups  were  detected  at  the  lowest  frequencies  (3%  and  9%, 
respectively;  Table  4). 

In  Honduras,  as  in  South  Africa,  the  Dcfe\  variant  group  with  the  highest 
frequency  (64%)  in  drones  was  D400  (Table  5).  D403,  D301 ,  and  D41 4,  were 
detected  at  the  highest  individual  frequencies  (19%,  13%,  and  13%, 
respectively).    Frequencies  were  low  for  all  other  variants  detected  (Table  5). 

Individual  frequencies  were  low  for  the  many  alleles  detected  in  drones 
from  Honduras,  as  in  the  South  African  drones  in  which  similar  alleles  were 
detected  (Table  6).  The  allele  detected  at  the  highest  frequency  in  Honduras 
(15%)  was  M301/D301,  the  allele  detected  at  the  highest  frequency  in  Mexico. 

Costa  Rica.  Examples  of  the  Msp\  restriction  fragment  patterns  detected 
in  the  Costa  Rican  colonies  maintained  at  2200m  are  shown  in  Figure  9,  lanes 
1,  2,  4-9,  and  12  (possibly  lane  14  as  well).  In  these  colonies  (Table  4,  Costa 
Rica,  'European'),  the  group  of  Msp\  variants  with  the  highest  frequency  was 
Ml 00  (38%),  followed  by  M300  (35%).  The  detection  of  M400  and  M500 
variants  indicated  that  African  paternal  introgression  had  occurred  in  these 
colonies.  The  M500  group  frequency  may  have  been  as  low  as  2%,  but 
fragment  comigration  confounded  the  differentiation  of  the  M200  and  M500 
variants  and  the  determination  of  their  frequencies.       ,_ 

Examples  of  the  restriction  fragment  patterns  detected  in  the  three  Costa 
Rican  colonies  maintained  at  700m  are  shown  in  Figure  9,  lanes  12,  14,  16, 
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and  17  (possibly  lane  9  as  well).  In  these  colonies  (Table  4,  Costa  Rica 
'African')  the  Msp\  variant  group  MSOO  was  detected  at  the  highest  frequency, 
54-58%.  The  M400  and  M300  group  frequencies  were  similar  to  each  other 
(Table  4).  . 

.  Venezuela.  Variants  from  all  five  Mspl  groups  were  detected  in  workers 

collected  in  Venezuela  (Figure  9,  lanes  4-6,  8,  9,  12,  14-17,  and  possibly  1,  3, 
and  11).  As  in  the  South  African  samples,  the  variant  group  detected  at  the 
highest  frequency  was  M500,  followed  by  the  M300  and  M400  groups.  M 1 00 
and  M200  variants  were  detected  at  the  lowest  frequencies  (Table  4).  The 
M400/M500  genotype  (Figure  9,  lane  16)  was  detected  most  frequently, 
followed  by  M500/M500  and  M300/M500  genotypes;  the  other  genotypes 
were  detected  at  much  lower  frequencies. 
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Old  World  bees.  v""'.     .'      >      .      ■4.^-.- 

In  the  east  European  workers  examined,  Msp\  variants  from  the  Ml 00 

group  appeared  to  be  specific  to  the  Italian  subspecies  ^4.  m.  ligustica,  while  in 

A.  m.  causcasica,  the  Msp\  restriction  fragment  patterns  consisted  only  of 

variants  belonging  to  the  M200  group.    Msp\  variants  from  the  M300  group 

were  exclusive  to  the  west  European  black  or  German  bee,  A.  m.  melHfera. 

The  Dde\  restriction  fragment  patterns  in  A.  m.  ligustica  and  A.  m.  caucasica 
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appeared  to  consist  only  of  variants  from  the  D300  group.    A.  m.  mellifera 

samples  examined  contained  Dde\  variants  from  the  D100,  D200,  and  D300 

groups.    These  distributions  are  not  absolute,  but  accurately  reflect  current 

sample  availability.        ► 

A  detailed  accounting  of  the  variation  in  A.  m.  scutellata  was  obtained 

from  the  identification  of  alleles  composed  of  Msp\  and  Dde\  variants  in  drones 

from  South  Africa,  collected  in  the  area  from  which  the  New  World  African 

bees  originated.  ;    ^      ■ 

North  American  bees. 

The  genetic  variation  at  locus  178  in  New  World  European  bees  was 
ascertained  by  examining  drone  bees  in  the  USA  and  feral  workers  from 
northern  Mexico.  The  Msp\  variant  groups  were  the  same  as  those  found  in 
Old  World  European  bees.  The  frequencies  for  the  Ml 00  restriction  fragment 
patterns  in  A.  m.  ligustica  workers  from  Europe  and  workers  from  the  USA 
were  nearly  identical.  In  USA  populations,  the  high  frequency  for  the  Ml 00 
variant  group  reflects  the  preferred  use  of  A.  m.  ligustica  for  beekeeping. 
These  results  concur  with  the  frequencies  in  USA  bees  reported  for  east 
European-specific  nuclear  and  mitochondrial  DNA  markers  (Hall  1986,  1990; 
Hall  &  Smith  1991;  Schiff  &  Sheppard  1993),  as  well  as  for  malate 
dehydrogenase  (MDH)  allozymes  (Nunamaker,  Wilson  &  Haley  1984;  Sheppard 
1988).   The  low  frequency  (9%)  for  the  Ml 00  variant  group  in  bees  collected 
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in  northern  Mexico  may  be  a  consequence  of  small  sample  size,  but  more  likely 

relects  the  small  east  European  contribution  to  the  feral  population  in  this  area. 

;y  Frequencies  for  the  Msp\  variant  group  M300  in  the  USA  (18%)  and  in 

northern  Mexico,  prior  to  the  arrival  of  African  bees  (64%  -  the  highest  of  all 

the  New  World  populations  examined),  are  consistent  with  the  history  of  the 

importation  of  west  European  bees  in  North  America,  their  limited  commercial 

use,  and  persistence  as  feral  colonies.    The  presence  of  the  west  European 

variant  M301  is  also  in  agreement  with  MDH  allozyme  frequencies  in  USA  bees 

(which  indicate  >4.  m.  mellifera  nuclear  genes  persist  in  feral  colonies;  Sheppard 

1 988,  1 989).  West  European  mitochondrial  DNA  was  detected  at  frequencies 

of  7%  in  USA  managed  colonies  and  64%  in  feral  colonies  in  northern  Mexico 

prior  to  the  arrival  of  African  bees  (Hall  &  Smith  1991)  and  was  detected  in 

21  %  of  422  feral  colonies  in  the  south-central  and  southeastern  USA  (Schiff 

&  Sheppard   1993).     Thus,   bees  of  west  European  ancestry  continue  to 

contribute  to  the  honey  bee  gene  pool  in  feral  colonies  of  North  America. 

Dde\  variant  group  D300  was  detected  at  the  highest  frequency  in  USA 

bees,  consistent  with  the  detection  of  variants  in  this  group  in  each  of  the  three 

European  races  examined,  and  at  the  highest  frequency  in    A.  m.  ligustica. 

Dde\  variants  D401  and  D404  were  found  in  two  separate  USA  colonies;  D401 

was  found  in  a  feral  colony  near  Tucson,  Arizona,  and  D404  was  found  in  a 

managed  colony  in  Kansas.    This  finding  was  noteworthy  because  the  D400 

variant  group  was  found  at  high  frequency  in  South  African  bees.    D401  was 
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not  found  in  the  South  African  colonies,  and  D404  was  detected  in  a  single 

South  African  drone.  One  explanation  for  the  detection  of  these  variants  in  the 

USA  is  that  the  distribution  of  the  Dde\  variants  may  not  be  absolutely 

discontinuous  and  correlated  with  fragment  pattern.     Alternatively,  these 

variants  may  indicate  African  ancestry  at  low  frequency  in  USA  bees.  Evidence 

for  the  persistence  of  A.  m.  lamarkii  in  the  USA,  which  was  introduced  (and 

soon  abandoned)  for  beekeeping  in  1866  (Schiff  &  Sheppard  1993;  Sheppard 

1989)  was  recently  shown  by  the  detection  of  non-scutellata  African  mtDNA 

at  a  frequency  of  1  %  in  feral  bees  in  the  south-central  and  southeastern  USA 

(Schiff  &  Sheppard  1993).  Variants  D401  and  D404  may  reflect /I.  m.  lamarkii 

ancestry,  but  this  possibility  needs  to  be  confirmed  with  samples  of  this  race. 

On  the  basis  of  the  congruity  in  the  Mspl  and  Dcle\  variants  detected  in 

Europe  and  in  North  America,  the  alleles  identified  in  drones  in  the  USA 

represented  a  subset  of  the  variability  in  European  bees.  The  allele  detected  at 

the  highest  frequency  (39%)  in  USA  drones,  M 1 03/D303,  may  be  specific  for 

A.  m.  ligustica  ancestry.  In  USA  bees,  only  variants  from  the  D300  group  were 

detected   with   Ml 00   and   M200   variants,   consistent  with   the  restriction 

fragment  patterns  in  the  east  European  samples.    Mspl  variants  M301   and 

M302  in  USA  bees  were  detected  with  Dde\  variant  D301  and  variants  in  the 

D1 00  and  D200  groups,  evidence  of  which  was  also  seen  in  the  A.  m.  mellifera 

workers  in  which  only  M300  variants  were  present.    Alleles  M101/D401  and 
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M302/D404  were  present  at  very  low  frequency;  these  Dde\  variants,  but  not 

these  alleles,  were  also  found  in  South  African  bees. 
Neotropical  bees. 

■"'•*,. 

'  ■'*)■■■ 

Italian  and  German  farmers  brought /I.  m.  ligustica  and  A.  m.  meff/fera, 
respectively,  to  their  settlements  in  south  and  southeastern  Brazil  (Gongalves 
1 974;  Goncalves,  Stort  &  DeJong  1 991 ;  Lobo,  Del  Lama  &  Mestiner  1 989)  and 
Argentina  (Kerr,  De  Leon  &  Dardo  1982).  A.  m.  caucasica  and  >4.  m.  carnica 
were  also  introduced  (Kent  1988;  Ruttner  1986).  In  contrast  to  the  USA,  the 
contribution  of  A.  m.  ligustica  to  the  total  gene  pool  in  the  neotropics  was 
minor  in  spite  of  considerable  importation  (Kent  1988).  The  Italian  bee  has 
been  found  concentrated  in  certain  areas,  particularly  Argentina  (Dietz,  Krell  & 
Eischen  1985;  Kerr,  De  Leon  &  Dardo  1982;  Sheppard  ef  a/.  1991),  Costa  Rica 
(Hall  1990;  Kent  1988;  Spivak  1991),  and  the  Yucatan  (Kent  1988;  Rinderer 
et  al.  1991).  There  is  little  indication  that  significant  introductions  of  east 
European  bees  were  made  or  were  successful  in  central  and  northeastern  Brazil 
(Lobo,  Del  Lama  &  Mestriner  1989),  the  Guianas,  Suriname  (Taylor  1977),  or 
Panama  (Boreham  &  Roubik  1 987;  Roubik  1 982).  Prior  to  the  release  of  A.  m. 
scute/fata   in   South   America,   A.    m.    mellifera   remained   dominant  in   the 


-^  neotropics,  and  it  has  been  assumed  that  if  early  introductions  of  subspecies 

originating  outside  Europe  were  made,  they  were  unsuccessful  (Kent  1988). 
Frequencies  for  the   variants   and   variant  groups  identified   in   east 
European  workers  and  in  the  USA  were  very  low,  or  absent,  in  neotropical  bees 
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collected  in  areas  where  African  bees  were  established.    The  east  European 

variants  were  largely  in  managed  colonies.  In  southern  Mexico,  the  frequency 
for  the  Ml 00  group  of  variants  was  24%  in  drones,  and  10%  in  workers. 
These  frequencies  are  consistent  with  previously  identified  nuclear  DNA  RFLP 
markers  specific  to  east  European  bees  (Hall  1986,  1990),  and  may  have 
reflected  more  recent  colonization  by  African  bees  compared  to  Honduras  and 
Venezuela.  Some  of  these  managed  colonies  were  known  to  be  of  east 
European  ancestry,  based  on  mtDNA  (Hall  &  Muralidharan  1989;  Hall  &  Smith 
1991;  Smith,  Taylor  &  Brown  1989).  ,   ,      .^: 

The  detection  of  the  Msp\  variant  M301  in  colonies  established  primarily 
from  feral  swarms  provided  evidence  for  the  persistence  of  A.  m.  mellifera  in 
the  neotropics,  as  has  been  suggested  by  allozyme  frequencies  (Lobo,  Del  Lama 
&  Mestriner  1989).  The  M300  group  frequency  ranged  from  19-33%, 
including  areas  which  have  been  occupied  by  African  bees  for  over  20  years. 

In  workers  from  Honduras,  the  detection  of  an  Msp\  variant  bearing  the 
1 .8kb  fragment  could  not  be  identified  outright  as  M301  and  thus  attributed  to 
A.  m.  mellifera  ancestry  in  the  apiaries  and  feral  populations  from  which  the 
samples  were  obtained,  since  variants  M303  and  M304,  found  in  South  African 
drones  at  much  lower  frequencies,  also  contained  this  fragment  (Figure  8; 
Tables  3,  4  &  6).  The  possibility  of  African  origin  for  the  1 .8kb  Msp\  fragment 
had  to  be  considered,  because  other  Msp\  variants  had  been  detected  at  higher 
frequencies  in  Honduras  than  in  South  Africa  (e.g.,  M400  variants).  It  was  also 

-  -    .'■  V 


64 
possible  that  a  variant  bearing  the  1.8kb  Msp\  fragment  detected  in  South 

African  bees  may  have  originated  from  prior  Importation  of  bees  to  South  Africa 

from  Europe  or  by  way  of  North  America  (Fletcher  1973,  1978). 

The  west  European  origin  of  the  I.Skb  fragment-bearing  Msp\  variants 

(M300   group)   in   southern   Mexico   and   Honduras  was   confirmed   by  the 

identification  of  the  associated  Dde\  variants  at  locus  1 78  (Table  6,  Figure  11). 

In  South  African  drones  the  Msp\  M303  and  M304  variants  were  associated 

with  the  Dde\  variants  D405  and  D503,  respectively.    In  neotropical  drones, 

M301   was  identified  as  the  Msp\  variant  with  the  I.Skb  fragment  and  was 

associated  with  Dde\  variants  that  were  the  same  as  (D101,  D102,  D201, 

D301 )  or  similar  to  (D202,  D203,  D305)  Dde\  variants  identified  in  USA  drones 

and  inferred  from  Dde\  fragments  detected  in  A.  m.  mellifera  workers.  Among 

the  M300  group  variants,  the  detection  of  only  M301  and  M302  in  Honduras 

demonstrates  that  A.  m.  mellifera  ancestry  persists  in  feral  colonies.    Dde\ 

D200  variants  detected  with  M301  in  Honduras  may  indicate  A.  m.  mellifera 

ancestry,  but  fragment  comigration  prohibited  the  detection  of  these  Dde\ 

variants  in  A.  m.  mellifera  workers.   There  was  no  evidence  for  Dde\  variants 

D405  or  D503  in  A.  m.  mellifera  or  in  the  M300-containing  samples  in  the 

neotropics.   The  absence  of  variants  M303  and  M304  In  the  New  World  may 

be  due  to  sampling  error;  these  variants  may  not  have  been  present  in  the  A. 

m.  scutellata  queens  imported,  or  were  present  at  a  low  frequency  not  detected 

in  this  study.  ■^.  ' 
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LOCALIZATION  OF  VARIATION  AT  LOCUS  178  IN  Apis  mellifera  (L.) 

BY  RESTRICTION  MAPPING 


Introduction 

African  and  European  groups  of  honey  bee  subspecies  have  been  shown 
to  differ  at  locus  178  in  the  number  and  location  of  the  four-  and  five- 
nucleotide  base  recognition  sequences  of  the  restriction  enzymes  Msp\  and 
Dde\,  respectively.  While  the  ultimate  resolution  of  the  variation  at  locus  178 
occurs  at  the  level  of  the  DNA  sequence,  RFLP  analysis  has  permitted  an 
assessment  of  the  variation  using  a  small  portion  of  its  nucleotides. 

The  five  groups  of  variants  identified  at  178  for  both  Msp\  and  Dcle\, 
which  have  similar  restriction  fragments  and  geographic  distributions,  may 
share  a  subset  of  restriction  sites.  The  number  of  fragments  present  in  all 
variants  within  a  group  ranged  from  one  in  the  M500  group  to  ten  in  the  M200 
group  (Figure  8).  Between  any  two  groups  of  Msp\  or  Dde\  variants,  the 
number  of  common  fragments  ranged  from  one  to  seven.  One  fragment  was 
common  to  all  Msp\  variants  (Figure  8),  and  no  fragments  appeared  to  be 
common  to  all  Dcle\  variants  (Figure  11).  It  is  likely  that  a  specific  composite 
of  restriction  sites  exits  for  each  subspecies  or  group  of  subspecies. 
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RFLP  analysis  has  revealed  greater  subspecies-specific  genetic  variation 

in  honey  bees  than  has  any  other  method  to  date.  However,  for  routine  use 
(e.g.,  identification  for  regulation,  large  scale  population  studies),  a  more  rapid 
method  of  analysis  would  be  required.  Limitations  of  RFLP  analysis  for  routine 
use  include  the  isolation  of  sufficient  genomic  DNA  for  two  restriction 
digests/Southern  blots,  the  (almost  certain)  use  of  radioisotope  for  labeling  the 
probe  DNA,  and  the  amount  of  time  for  obtaining  results.  In  addition,  a 
different  approach  was  needed  to  investigate  the  allelic  nature  of  the  restriction 
fragments  in  each  digest,  and  to  identify  subspecies-specific  length 
polymorphisms  and/or  polymorphic  restriction  sites.  For  such  objectives,  the 
slight  differences  between  variants  in  the  sizes  of  restriction  fragments,  the 
complex  nature  of  the  restriction  fragment  patterns,  and  the  size  of  the  locus 
precluded  the  continued  use  of  the  RFLP/Southern  blot  technique. 

Conversion  of  the  RFLP  analysis  to  a  PCR  (polymerase  chain  reaction)- 
compatible  format  (Saiki  et  a/.,  1988)  would  expedite  the  identification  of 
individual  bees,  the  determination  of  the  allelic  nature  of  the  restriction 
fragments,  and  permit  the  use  of  the  polymorphic  sites  to  initiate  investigations 
of  the  variation,  organization,  and  genetic  content  of  locus  178.  The 
localization  of  restriction  site  and  length  polymorphisms  to  specific  regions  of 
the  locus  would  identify  the  region(s)  for  which  PCR-primers  could  be  made. 
The  conversion  process  was  initiated  by  isolating  smaller  regions  of  probe 
pB178,  identifying  the  Msp\  and  Dde\  restriction  fragments  to  which  each 
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region  hybridized  in  Southern  blots,  and  identifying  Msp\  and  Dde\  sites  within 

each  region. 


Materials  and  Methods 

The  cloned  probe  pBI  78  was  replicated  in  Escherichia  coli  strain  DH5ar 
(Gibco  BRL)  and  isolated  by  the  alkaline-lysis  procedure  (Sambrook,  Fritsch  & 
Maniatis  1989). 

Mapping  restriction  sites  within  pB178.  Purified  pB178  was  digested 
with  Pst\  to  release  the  9.45kb  honey  bee  DNA  insert  (referred  to  as  178), 
which  was  subsequently  separated  from  the  vector  by  electrophoresis  in  0.8% 
SeaPlaque  agarose  (FMC)  and  isolated  by  extraction  with  phenol  and  ether. 
The  relative  positions  of  the  restriction  sites  were  determined  by  sequentially 
and  reciprocally  digesting  178  and  pB178  with  up  to  four  restriction  enzymes. 
pB178  was  treated  with  one  or  more  restriction  enzymes,  and  one  aliquot  was 
set  aside  for  analysis  while  another  was  digested  with  Pst\  to  determine  the 
location  of  the  endonuclease  recognition  sites  relative  to  the  Pst\  ends  of  the 
insert.  The  digested  DNA  was  electrophoresed  in  1  %  agarose  (Kodak  IBI),  and 
the  fragments  were  visualized  by  staining  with  ethidium  bromide.  Fragment 
sizes  were  estimated  by  comparison  to  molecular  weight  standards  using  a  HI- 
PAD  digitizer.  For  comparison  and  confirmation,  the  fragment  sizes  were  also 
estimated  using  the  National  Center  for  Supercomputing  Applications  (NCSA) 
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GelReader  2.0  shareware  package  for  the  Macintosh  (NCSA  Software  Tool 

Group  at  the  Center  for  Prokaryotic  Genome  Analysis,  University  of  Illinois). 

Mapping  Msp\  and  Dde\  sites.  Smaller  regions  of  1 78  were  purified  from 
SeaPlaque.  These  regions  were  digested  with  Msp\  and  Dde\  in  the  presence 
and  absence  of  other  restriction  enzymes  known  to  have  recognition  sites 
within  the  region.  One-half  of  each  digest  was  electrophoresed  in  3%  agarose 
and  visualized  by  staining  with  ethidium  bromide;  the  other  half  was  end- 
labeled  (Sambrook,  Fritsch  &  Maniatis  1989),  electrophoresed  in  agarose,  and 
visualized  by  autoradiography.  Fragment  sizes  were  estimated  using  the 
GelReader  program.  ,   -      '     ' 

The  correspondence  of  the  fragments  composing  the  Msp\  and  Dde\ 
variants  to  the  physical  map  was  determined  by  probing  Southern  blots  with 
smaller  regions  of  178.  The  same  regions  of  178  used  for  mapping  the  Msp\ 
and  Dde\  sites  were  sequentially  radiolabeled  by  random-priming  (Feinberg  & 
Vogelstein  1983,  1984)  and  hybridized  to  the  blots  shown  in  Figures  2,  3,  and 
5-7. 

Subcloning  pB178.  Most  of  the  9.45kb  insert  was  subcloned  into  the 
plasmid  vector  pGEM3Z  (Promega).  The  subclones  were  given  the  prefix 
pG178,  and  each  was  identified  by  a  suffix  specific  for  the  region  subcloned 
(see  Figure  1 2).  The  pGI  78  subclones  were  replicated  and  purified  in  the  same 
manner  as  pB178.  The  subcloned  inserts  were  excised  and  extracted  from 
SeaPlaque.  The  ends  of  the  subclones  were  sequenced,  using  the  M 1 3  forward 
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;  and  reverse  sequencing  primer  sites  in  the  vector,  at  the  University  of  Florida 

Interdisciplinary  Center  for  Biotechnology  Research  sequencing  facility.    The 
sequences  were  checked  for  Msp\  and  Dde\  sites  using  the  Seqaid  II  program 
r        version  3.5  (D.  Rhoads  &  D.  Roufa,  Kansas  State  University). 

r;  .,  '  Results 

f».'  Restriction  sites  in  178  used  for  subcloning,  for  mapping  Mspl  and  Df/el 

sites,  and  for  Southern  blot  hybridizations  are  shown  in  Figure  12.    Five  non- 

-  y  '"T 

'    i      overlapping  fragments,  representing  the  entire  178  insert,  were  isolated  and 
^,:  used  to  probe  Southern  blots  (Figure  12):    178PiK,  178KH,  178HB,  I78BE4, 

and  I78E4P2.  The  latter  four  regions  were  subcloned  (Figure  12)  and  the 
inserts  were  used  to  map  Msp\  and  Dde\  sites.  The  2.6kb  Pst\^  to  Kpn\ 
fragment  (1 78PiK)  was  isolated  from  pBI  78  but  was  not  subcloned.  The  first 
l.lkbof  178PiK,  fromPsfli  toEcoRIi  (178PiEi,  Figure  12),  was  subcloned  and 
the  Msp\  and  Dde\  sites  were  determined,  but  it  was  not  used  to  probe 
Southern  blots. 
Identification  of  Msp\  sites  in  178.  *  ■  i     -l 

The  Msp\   sites  identified  in    1 78   are  shown  in   Figure   1 3  and  are 
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N  presented  in  the  same  format  as  Figure  12.    The  results  of  the  Southern  blot 

hybridizations,  shown  in  Figures  14  -  19,  are  presented  in  the  same  format  as 
Figure  1.      . 
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Figure  1 2.  Map  of  restriction  sites  identified  within  the  9.45kb,  fsfl-honey  bee 
genomic  DNA  insert  of  probe  pB178.  Nunnbers  to  the  left  of  the  continuous 
verticle  line  of  the  map  indicate  the  estimated  position  of  the  restriction  site 
(identified  on  the  right)  relative  to  one  of  the  Pst\  insertion  sites  (PsflJ. 
Relative  positions  of  sites  were  determined  by  digesting  the  intact  clone  pBI  78 
with  each  enzyme  in  the  presence  and  in  the  absence  of  Pst\,  and  in 
combination  with  up  to  three  additional  enzymes.  These  sites  were  confirmed 
by  isolating  the  9.45kb  insert  (digestion  of  pB178  with  Pst\  followed  by 
electrophoresis  in  and  extraction  from  0.8%  Seaplaque)  and  repeating  the 
single  and  multiple  enzyme  digests  in  the  absence  of  Pst\.  These  sites  served 
as  landmarks  for  determining  the  locations  of  the  Msp\  sites  (Figure  1 3)  and  for 
the  Dt/el  sites  (Figure  20)  within  the  probe. 

The  short  verticle  lines  on  the  left  side  of  the  page,  together  with  the 
letters  in  bold,  indicate  the  regions  of  pB178  subcloned  into  pGEM3Z.  The 
pGEMSZ  subclones  were  identified  according  to  the  region  contained  in  each: 

pG178PiEi:    l.lkb,  Psfl,  site  (top)  to  the  first  £coRI  site; 

pG178KH:     0.9kb, /CpA7l  -  M>7dlll; 

PG178HB:      :  .9kb,  Hind\\\  -  Bg/W; 

PGI78BE4:     2.0kb,  Bp/ll  -  EcoRI; 

PGI78S3E4:  0.75kb,  third  Sa/I  site  to  the  fourth  EcoRI  site; 

PGI78E4P2:  2.0kb,  fourth  £coRI  site  to  terminal  Psflj  site. 

The  dotted  verticle  line  indicates  the  2.6kb  Pst\^-Kpn\  region  which  was 
used  for  mapping  and  hybridized  to  Southern  blots,  but  was  not  subcloned. 
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Digestion  of  178PiK  with  Msp\  produced  fragments  of  approximately 

0.73kb,  0.62l<b,  0.41  kb,  and  0.38kb  (Figure  1 3).  Double  digestion  of  1 78PiK 
with  Mspl  and  each  restriction  enzyme  with  a  recognition  site  within  178PiK 
(Figure  12)  revealed  that  there  were  two  0.38kb  Mspl  fragments.  The 
fragments  to  which  178PiK  hybridized  are  shown  in  Figure  14.  The  intense 
hybridization  signal  seen  at0.38kb  in  the  Ml 00  variants  in  the  hybridization  of 
both  pBI  78  and  1 78PiK  was  due  to  the  presence  of  two  fragments.  Most  of 
the  variants  were  identical  in  the  fragments  to  which  178PiK  hybridized.  The 
Mspl  restriction  fragment  pattern  of  178PiK  resembled  Mspl  variants  M101 
and  Ml  03.  Additional  Mspl  sites,  shown  in  Figure  1 5,  have  been  hypothesized 
to  account  for  the  fragments  seen  in  other  variants. 

Digestion  of  the  0.9kb  Kpn\  -  HindWl  region  of  178  (178KH)  with  Mspl 
produced  a  single  visible  fragment  of  approximately  0.65kb.  Two  Mspl  sites 
were  identified  near  the  Kpn\  site  in  the  sequence  of  the  pG178KH  subclone 
(Figure  13:  these  fragments  had  not  been  detected  following  the  digests  due 
to  the  small  amount  of  DNA  used  in  the  gels).  The  Mspl  fragments  to  which 
1 78KH  hybridized  are  shown  in  Figure  1 6.  In  the  majority  of  variants,  1 78KH 
hybridized  to  a  0.75kb  fragment  or  to  one  just  slightly  larger  or  smaller.  In 
these  variants  a  very  faint  hybridization  signal  was  seen  at  0.56kb  (not  shown 
in  Figure  16),  which  was  also  detected  by  178PiK.  In  variants  Ml 02  and 
M201,  178KH  hybridized  to  a  0.75kb  fragment,  as  well  as  to  a  0.62kb 
fragment  to  which  1 78P,K  appeared  to  have  hybridized.  As  indicated  in  Figure 


Figure  13.  Map  of  Msp\  restriction  sites  identified  in  honey  bee  genomic  DNA 
insert  of  probe  pB178.  The  Pst\  sites  marking  the  ends  of  the  insert  were 
utilized  for  cloning  in  pBR322.  Relative  positions  of  sites  were  determined  by 
digesting  smaller  regions  of  the  insert  (Pst\^  -  Kpn\,  Kpn\  -  H/nd\\\,  Hind\\\  -  Bg/\\, 
Bg/\\  -  £coRl4,  and  EcoR\^  -Psflj)  with  Msp\  in  the  absence  and  presence  of 
restriction  enzymes  within  each  region,  which  were  indicated  in  Figure  12 
(with  the  exception  of  the  region  from  the  Kpn\  site  to  the  H/nd\\\  site).  The 
asterisks  (*)  indicate  sites  identified  in  the  DNA  sequences  of  the  ends  of  the 
respective  subclones.   Left:   distances  (in  base  pairs,  bp)  between  Msp\  sites. 
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:{';^;  ,        16,  the  loss  of  an  Msp\  site  in  this  region  (at  approximately  2785bp)  could 

explain  the  presence  of  the  0.62kb  fragment  in  these  two  variants. 

Digestion  of  the  1.9kb  region  from  Hin6\\\  to  Bgl\\  (178HB)  with  Msp\ 
produced  fragments  of  approximately  1.3kb,  0.42kb,  lOObp  and  70bp.   Msp\ 

^,.  sites  were  identified  near  the  Hin6\\\  and  Bgl\\  sites  in  subclone  pG178HB 

*  (Figure  1 3).   Msp\  restriction  fragments  to  which  1 78HB  hybridized  are  shown 

in  Figure  1 7.  The  Msp\  restriction  fragment  pattern  of  1 78HB  corresponded  to 
variant  Ml 03:    two  site  polymorphisms  are  indicated  in  Figure  17  that  could 
; '  account  for  the  fragments  detected  by  1 78HB  in  two  other  variants. 

I;;;  Digestion  of  the  2.0kb  Bgl\\  -  EcoRU  fragment  (I78BE4)  with  Msp\ 

produced  fragments  of  approximately  1.4kb,  0.33kb,  and  0.25kb  (Figure  13). 
The  Msp\  fragments  to  which  1 78BE4  hybridized  in  Southern  blots  are  shown 
in  Figure  18:  these  included  combinations  of  approximately  1.8kb,  1.4kb, 
1.25kb,  0.77kb,  0.56kb,  and  0.35kb  fragments.  The  loss  of  an  Msp\  site  in 
the  178HB  region  (Figures  13  &  17,  at  approximately  5420bp)  could  account 
for  the  0.77kb  fragment  detected  in  variants  Ml 01  and  M302,  and  for  the 
detection  of  this  fragment  with  both  178HB  and  I78BE4.  Since  I78BE4 
/  hybridized  to  more  Msp\  fragments  than  had  been  found  in  the  digestion  of  the 

insert  and  accounted  for  by  overlapping  with  adjacent  regions,  the  placement 
of  Msp\  sites  within  1 78BE4  was  based  on  double  digests  with  Spe\,  Sph\,  Nsi\, 

^  and  Sa/I,  and  the  results  of  the  hybridizations  of  178HB,  I78BE4,  and  I78E4P2 

V.^;'         (following). 
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Digestion  of  the  2.0kb  EcoR\^  -  Pst\2  fragment  (I78E4P2)  with  Msp\ 
produced  fragments  of  approximately  1 .3kb  and  0.56kb.  Two  additional  Msp\ 
sites  were  found  near  the  EcoR\  site  in  subclone  PGI78E4P2  (Figure  13).  The 
fragments  to  which  I78E4P2  hybridized,  shown  in  Figure  19,  appeared  to 
include  the  1  .Bkb,  1 .25kb,  and  0.56kb  fragments  to  which  1 78BE4  hybridized. 
In  the  majority  of  variants,  178E4P2  hybridized  to  fragments  of  approximately 
1 .25kb  and  0.56kb.  The  existence  of  the  1 .8kb  fragment  in  the  M300  variants 
could  be  explained  by  the  loss  of  the  Msp\  site  at  8140bp  (Figure  13).         ''  — 

Since  it  is  likely  that  the  duplication  in  the  fragments  to  which  1 78BE4 
and  1 78E4P2  hybridized  resulted  from  contamination  of  1 78BE4  with  1 78E4P2, 
the  placement  of  the  Msp\  sites  shown  in  Figure  1 3  was  based  on  the  digestion 
of  different  aliquots  of  these  two  regions  with  Msp\,  in  the  absence  and 
presence  of  other  restriction  enzymes  which  have  sites  in  each  region. 

Identification  of  Dde\  sites  in  1 78.  ■  , 

The  approximate  locations  of  Dde\  sites  identified  in  the  honey  bee  DNA 
insert  of  pBI  78  are  shown  in  Figure  20,  in  the  same  format  as  Figures  1 2  and 
13.  The  results  of  Southern  blot  hybridizations  are  shown  in  Figures  21  -  25, 
in  the  same  format  as  Figure  4.  ^  ; '■ 

Digestion  of  178PiK  with  Dcle\  produced  fragments  of  approximately 
1.3kb,  0.75kb,  and  0.4kb.  A  single  fragment  of  approximately  I.Okb  was 
detected  when  178PiEi  was  digested  with  Dde\  (Figure  20:   refer  to  Figure  13 


A 


,i  '■"■  t 


Figure  20.  Map  of  Dde\  restriction  sites  identified  in  honey  bee  genomic  DNA 
insert  of  probe  pB178.  The  Pst\  sites  marl<ing  the  ends  of  the  insert  were 
utilized  for  cloning  in  pBR322.  Relative  positions  of  sites  were  deternnined  by 
digesting  snnaller  regions  of  the  insert  (Pst\^  -  Kpn\,  Kpn\  -  HindWl,  HindlW  -  Bg/\\, 
Bg/\\  -  EcoRl^,  and  EcoR\^  -  Psflj)  with  Dde\  in  the  absence  and  presence  of 
restriction  enzymes  within  each  region,  which  were  indicated  in  Figure  1 2  (with 
the  exception  of  the  region  between  the  Kpn\  and  HindWl  sites).  The  asterisks 
(*)  indicate  sites  identified  by  sequencing  the  ends  of  the  respective  subclones. 
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V      86 
for  178PiEi  location  and  size).     Analysis  of  the  sequences  of  the  ends  of 

subclone  pG178PiEi  revealed  a  Dde\  site  at  150bp  from  the  Pst\^  site  (Figures 

20  &  21 ).  The  Dde\  fragnnents  to  which  1 78PiK  hybridized  are  shown  in  Figure 

21.       In   the    nnajority    of   variants,    178PiK   hybridized   to   a   fragment   of 

;,  approximately  0.75kb,  and  another  fragment  of  approximately  I.Skb.    In  the 

%■'/         D100  and  D200  variants,  a  single  fragment  of  approximately  2.1kb  was 

detected.   The  loss  of  the  Dde\  site  at  approximately  1450bp  would  account 

for  the  2.1kb  fragment  in  the  D100  and  D200  variants.    The  variation  in  the 

size  of  the  fragment  of  approximately  1.3kb  may  be  the  result  of  length 

polymorphism(s)  closer  to  the  Psfli  site.  -     '^ 

; '"  :  -  ■\':^-^--  -j^f-v-  -'  V' 

There  were  no  Dde\  sites  within  178KH.     A  single  fragment  in  each 

■>  ■'■ 

>W  variant,  approximately  2.0kb,  was  detected  when  178KH  was  used  to  probe 
this  region,  shown  in  Figure  22.  On  the  basis  of  the  Dde\  sites  identified  in 
1 78PiK,  it  was  concluded  that  the  same  2.0kb  fragment  had  been  detected  by 
1 78PiK  and  1 78KH,  and  would  be  detected  by  1 78HB  as  well.  The  placement 
of  the  Dde\  sites  on  either  side  of  178KH  is  shown  in  Figure  22. 

Digestion  of  1 78HB  with  Dde\  resulted  in  fragments  of  approximately 
1.1 5kb  and  0.75kb.  The  Dde\  fragments  to  which  178HB  hybridized  are 
shown  in  Figure  23.  In  each  variant,  178HB  hybridized  to  the  2.0kb  fragment 
detected  by  1 78PiK  and  1 78KH.  One  other  fragment  was  detected  by  1 78HB 
in  each  variant:  approximately  I.Okb  in  the  DIOOs  and  D400s,  1.25kb  in  the 
D200s  and  D300s,  and  0.98kb  in  the  D500s.  Two  possible  locations  for  a  site 
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^1  ■        which  could  account  for  the  I.Okb  fragment  seen  in  the  D100  and  D400 
»||  variants  are  shown  in  Figure  23.  /*;... 

Digestion  of  178BE4  with  Dde\  resulted  in  fragments  of  approximately 

1 .27kb,  0.39kb,  and  0.25kb.  The  Dde\  fragments  to  which  1 78BE4  hybridized 

'^       are  shown  in  Figure  24.     In  the  majority  of  variants  I78BE4  hybridized  to 

fragments  of  approximately  2.1kb,   I.Skb,  and  1.05kb.     A  Dde\  site  was 

■■>•■■■■ 

'•;  identified  near  the  Bgf\\  site  of  subclone  PGI78BE4  (Figure  20).   Since  I78BE4 

7k  hybridized  to  more  Dde\  fragments  than  it  appeared  to  contain,  which  also 

occurred  when  1 78BE4  was  used  to  probe  Southern  blots  containing  the  Msp\ 
variants,  the  placement  of  Msp\  sites  within  I78BE4  was  based  on  double 
-  digests  with  Spe\,  Sph\,  Nsi\,  and  Sa/I,  and  the  results  of  the  hybridizations  of 

178HB,  I78BE4,  and  I78E4P2  (following).  ^4       %:->4-^    .  '- 

Digestion  of  I78E4P2  with  Dde\  produced  fragments  of  approximately 
1 .2kb  and  0.85kb  (Figure  20).  The  fragments  to  which  1 78E4P2  hybridized  are    . 
shown  in  Figure  25,  and  include  what  appeared  to  be  the  same  2.1kb  and  " 
1 .05kb  fragments  to  which  1 78BE4  had  hybridized.  Sites  corresponding  to  the 
fragments  greater  than  2.3kb  in  the  D400  variants  lie  outside  the  probe  region. 

The  duplication  in  the  Dde\  fragments  to  which  I78BE4  and  I78E4P2   -:'■'■ 
hybridized  was  concluded  to  be  due  to  contamination  of  1 78BE4  with  1 78E4P2. 
Therefore  the  placement  of  the  Dde\  sites  shown  in  Figure  20  was  based  on  the 
digestion  of  different  aliquots  of  these  two  regions  with  Dde\,  in  the  absence 

V;  V.        and  presence  of  other  restriction  enzymes  with  sites  in  each  region. 

■'iV.  :    ■  ■■■■      -  :-,,v--        -  • "      ,^  - 

'J''^:  ■■■'•-.,. 
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The  sizes  of  fragments  detected  in  Southern  blots  and  by  ethidium 

bromide-staining  or  end-labeling  did  not  always  correspond.  The  migration  of 
DNA  molecules  in  agarose,  and  ultimately  the  determination  of  the  sizes  of  the 
DNA  fragments  so  separated,  depends  on  the  buffer  system,  the  voltage 
gradient,  and  the  concentration  of  the  agarose  (Southern  1979).  While  the 
same  buffer  was  used  for  all  electrophoretic  separations  reported  here,  the 
voltage  gradient  and  agarose  concentration  were  not  always  the  same.  All 
Southern  blots  involved  2%  agarose,  but  the  fragments  produced  from  the 
digestion  of  1 78  and  smaller  regions  of  1 78  were  separated  in  anywhere  from 
1  %  -  3%  agarose.  Although  more  than  135  Southern  blots  were  run,  some  of 

■■,■»■  "."■■",  ■ 

,■■:.-■  ■■  / 

the  fragment  sizes  from  the  Msp\  single  and  double  digestions  were  only 
determined  in  one  or  two  gels.  A  more  precise  estimation  of  the  sizes  of  the 
Msp\  and  Dde\  fragments  within  each  region  of  1 78  could  be  obtained  by 
repeating  the  digests,  gels,  and  measuring  the  fragment  migration  distances, 
followed  by  linear  regression  analysis.  However,  the  purpose  of  the 
hybridizations  and  digestions  performed  was  to  localize  polymorphisms,  which 
was  accomplished.  Further  details  can  be  obtained  from  the  proposed  PCR 
work  and  by  sequencing. 

;.-  Discussion 

The  Msp\  fragment  pattern  resulting  from  the  digestion  of  1 78  with  Msp\ 
resembled  the  pattern  of  Ml 03,  which  was  found  at  high  frequency  in  USA 


■,*fwj     ■    *._»-* 


bees  and  belonged  to  a  group  of  variants  found  at  very  high  frequency  in  bees 
sampled  from  east  Europe  and  North  America.  This  was  not  surprising,  as  the 
library  from  which  the  original  clone  pBI  78  was  obtained  was  constructed  with 
DN A  from  bees  of  east  European  ancestry  (Hall  1986). 

A  few  polymorphic  Msp\  sites  were  identified  in  178.  Additional  Msp\ 
sites  in  each  region  of  1 78,  not  detected  in  Msp\  digested  of  the  probe  DNA, 
have  been  hypothesized  to  account  for  the  fragments  detected  in   other 

:  ^vv  variants,  some  of  which  are  indicated  in  Figures  14  -  19.  For  example,  the 
M400  variants  as  a  group  were  characterized  by  a  unique  1 . 1  kb  fragment,  but 
lacked  a  0.85kb  fragment  present  in  all  of  the  other  Msp\  variants  (Figure  1). 
The  results  of  the  hybridization  with  178PiK  (Figure  15)  demonstrated  the 
allelic  relationship  of  the  1.1  kb  and  0.85kb  fragments. 

Ml 01,  Ml 03,  and  the  M400s  did  not  have  a  fragment  around  0.27kb 

which  was  observed  in  the  other  variants  (smallest  fragment  shown  in  Figure 

1).   The  Msp\  sites  which  could  account  for  the  1.1  kb  and  0.85kb  fragments 

-'-"■   ;       were  concluded  to  lie  outside  the  probe  region  (Figure  15,  sites  a  and  b),  in 

;:^:  which  case  the  Msp\  fragment  representing  the  difference  between  these  two 
sites  would  not  be  detected  by  the  probe.  The  site  that  accounts  for  the 
0.27kb  fragment  seen  in  all  variants  except  M101,  M103,  and  the  M400s 
could  then  be  located  between  sites  c  and  d,  or  between  sites  e  and  h.  The 
pattern  detected  for  M503  provided  the  answer;  it  contained  0.49kb  and 
0.27kb  fragments,  and  did  not  contain  a  0.38kb  fragment.    Site  f,  between 


U'.' 


m 

sites  e  and  h,  was  hypothesized  to  account  for  the  0.27kb  fragment;  the 
absence  of  site  d  in  M503  accounted  for  the  0.49kb  fragment;  and  an 
additional  site  g  was  hypothesized  to  account  for  the  absence  of  the  0.27kb 
fragment  in  the  M400  variants. 

In  Figures  1 ,  2,  and  3,  the  lower  molecular  weight  fragments  in  particular 
appear  to  increase  in  size  in  by  the  addition  of  a  sequence  of  constant  size. 
For  example,  some  of  the  additional  sites  hypothesized  in  regions  of  178 
(Figures  15,  16,  17,  &  19)  were  separated  by  approximately  55bp  or  llObp, 
or  some  combination  of  these  fragment  sizes.  178PiK  contained  two  O.SBkb 
fragments  (3  x  llObp  +  55bp)  separated  by  llObp,  as  well  as  a  0.56kb 
fragment  (5  x  1  lObp)  and  a  .62kb  fragment  (5  x  1  lObp  +  55bp).  The  loss  of 
a  site  in  the  1 78KH  region  accounts  for  the  allelic  relationship  between  the 
0.56kb  and  0.62kb  fragments  detected  in  different  variants  by  178PiK.  Msp\ 
sites  in  178PiK,  178KH,  and  178E4P2  resulted  in  three  different  locations  for 
0.56kb  fragments.  .  "^ 

A  polymorphic  Msp\  site  was  identified  in  I78E4P2,  at  approximately 
8140bp.  The  presence  of  this  site  was  manifest  in  Southern  blots  as  1.4kb 
and  0.56kb  fragments,  and  the  loss  of  the  site  resulted  in  a  1.8kb  fragment 
(Figures  13  &  19).  The  1 .8kb  fragment  was  characteristic  of  the  M300  group 
of  variants,  detected  at  a  frequency  of  94%  in  the  west  European  bee,  A.  m. 
mellifera;  at  frequencies  ranging  from  14%  -  64%  in  New  World  bees;  and  at 
a  frequency  of  3%  in  South  African  bees.     Individual  variants  in  the  group 
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appeared  to  be  geographically  restricted  in  distribution.  M302  was  detected 
in  bees  from  Europe  and  North  America,  and  was  distinguished  from  the  other 

■■^  three  variants  in  the  M300  group  by  the  presence  of  a  0.77kb  fragment,  which 

now  can  be  attributed  to  the  loss  of  an  Msp\  site  in  the  178HB  region. 

1^1     ?  M301   was  detected  in  bees  from  Europe,  USA,  and  the  neotropics; 

M303  and  M304  were  detected  only  in  South  Africa.  These  variants  could  be 

,  i'  difficult  to  distinguish  by  means  of  the  Msp\  fragment  patterns  alone,  but  were 

distinguished  by  the  Dde\  variants  associated  with  each  (Tables  2  &  6).  The 
ability  to  distinguish  the  alleles  containing  the  M300  variants  permited  the 
contribution  of  A.   m.   mellifera  to  the  feral  neotropical  populations  to  be 

:  '  ascertained,  which  had  not  been  possible  with  previously  found  nuclear  DNA 
markers.  The  regions  in  which  the  polymorphic  Msp\  and  Dde\  sites  are  located 
have  been  detemined,  although  the  sites  have  not  been  identified  on  a  fine 
scale.    Regions  have  been  identified  which  differentiate  between  most  of  the 

j>  Dde\  groups,  which  permits  a  preliminary  assessment,  presented  in  Table  7,  of 

:^;  •  the  regions  which  could  be  used  to  distinguish  the  alleles  containing  Msp\ 

M300  variants. 

c:    =  While  locus  178  represents  only  5  x  10"'%  of  the  honey  bee  genome 

( ~  1 .8x1  O^bp;  Jordan  &  Brosemer  1 974),  additional  analysis  of  locus  1 78  will 
provide  valuable  information  about  the  organization  of  the  genome,  about  which 
there  is  very  little:   only  mtDNA  has  been  studied  in  any  detail.   Such  analyses 

p,   •  will  initiate  a  data  base  for  determining  the  extent  of  variation  between 


i^ 


« 


97 


Table  7.    Polymorphic  Msp\  and  Dc/e\  restriction  sites  which  could  be  used  to 

distinguish  Msp\  M300  variants. 


M300  variant/ 
Dde\  variant  group 

D£yel-1450 

Restriction  site^ 
Mspl-5420            Dcrel-4530/5280 

M301/D200 

— 

+                                 — 

M301/D300 

+ 

+                                 — 

M302/D100              *: 

:  '  :— . '    ,"  -; 

;    "                 —                                                + 

M303/D400   ->.    ^r 
M304/D500 

+ 

■    ^    .^    *  ■'  .    U'i  ^^i?-^ •;>■'" 
+                                  + 

(  +  ),  Site  present;    ( — ),  site  absent 
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geographically  separated  populations  in  the  same  species,  and  markers  for 

following  the  movement  of  African  bees  in  the  USA. 

>    •         ,.       The  sequences  obtained  for  the  pGI  78  subclones  will  be  used  to  design 

primers  which  can  be  used  to  convert  the  RFLP  analysis  to  a  more  expedient 

PCR  format.     It  is  predicted  that  by  amplifying  each  region  in  the  samples 

containing    the    different    Msp\    and    Dde\    variants,    the    site    and    length 

polymophisms  which  define  each  variant  can  be  identified. 
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CONCLUDING  REMARKS 

RFLPs  at  the  polymorphic  locus  corresponding  to  probe  pBl  78  represent 
the  first  nuclear  DNA  markers  that  distinguish  east  European  or  north 
Mediterranean  (Garnery,  Cornuet  &  Solignac  1 992),  west  European,  and  South 
African  bees.  For  this  study,  only  workers  from  Europe  were  available,  and  the 
codominant  expression  and  comigration  of  restriction  fragments  in  diploid 
workers  precluded  the  identification  of  individual  variants  in  each  of  the 
European  populations.  Nonetheless,  the  detection  of  fragments  characteristic 
of  Msp\  and  Dde\  variant  groups  enabled  the  distribution  of  the  Msp\  and  Dde\ 
variants  in  the  European  samples  to  be  determined. 

A.  m.  mellifera  was  the  first  subspecies  introduced  to  North  America,  in 
1622  (Sheppard  1988),  followed  by /I.  m.  ligustica,  A.  m.  carnica,  and  A.  m. 
caucasica  in  the  mid-to-late  1800s  (Kent  1988;  Oertel  1976;  Pellet  1938; 
Sheppard  1989).  A.  m.  mellifera  overwintered  well  in  colder  climates,  and 
tended  to  establish  feral  colonies  upon  swarming,  but  this  race  was  defensive 
and  had  a  greater  susceptibility  to  brood  diseases  (Susceptibility  1982; 
Sheppard  1988).  The  east  European  bees,  particularly /4.  m.  ligustica,  gained 
tremendous  popularity  for  beekeeping  due  to  their  high  productivity  and  the 
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facility  of  their  management.    A.  m.  ligustica  has  dominated  the  hobby  and 
commercial  honey  bee  market  in  the  USA  for  over  a  hundred  years. 

It  is  likely  that  the  importation  of  honey  bees  from  Europe  to  the  USA 
resulted  in  a  reduction  in  variation,  or  genetic  bottleneck,  relative  to  the  parent 
populations  (Sheppard  1 988).  The  identity  and  number  of  variants  and  alleles 
at  locus  178  in  the  subspecies  imported  to  the  USA  for  beekeeping  have  not 
been  determined  but  may  be  greater  than  those  found  in  USA  drones.  Variants 
and  alleles  detected  in  USA  drones  likely  represent  a  subset  of  the  variability 
at  this  locus  in  European  honey  bee  subspecies.  The  emphasis  on  vitality  and 
productivity  in  beekeeping  and  queen-rearing  practices  has  resulted  in  some 
degree  of  homogenization  of  European  races  in  the  USA. 

Documentation  of  the  identity  and  numbers  of  European  honey  bee 
subspecies  imported  into  the  neotropics  appears  to  be  limited,  and  commercial 
and  feral  populations  have  not  been  thoroughly  surveyed  in  the  past  (Kent 
1988;  Taylor  1977).  It  has  been  reported  that  the  same  races  of  European 
bees  were  introduced  to  North  America  and  the  neotropics,  although  the 
introductions  occurred  at  different  times  (Goncalves  1974;  Goncalves,  Stort  & 
DeJong  1991;  Hellmich  &  Rinderer  1991;  Kent  1988;  Kerr,  DeLeon  &  Dardo 
1982;  Lobo,  Del  Lama  &  Mestriner  1989;  Rinderer  &  Hellmich  1991).  The 
number  of  African  queens  imported  and  released  in  Brazil  is  not  clear  (Kerr 
1967;  Rinderer,  Oldroyd  &  Sheppard  1993;  Smith  1991)  nor  the  number  of 
drones  with  which  each  was  mated,  but  the  introduction  probably  resulted  in 
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a  reduction  in  variation  relative  to  that  of  the  parent  South  African  population. 

Variants  and  alleles  found  in  the  South  African  samples  may  not  have  been 

represented  in  the  individuals  imported  to  the  New  World.    Variants  identified 

in  the  neotropical  samples  similar  to  but  not  detected  in  South  Africa  may 

reflect  sampling  error. 

Fletcher  has  provided  a  summary  of  past,  repeated  attempts  to  introduce 
A.  m.  ligustica  and  A.  m.  caucasica  to  South  Africa  (Fletcher  1973,  1978). 
The  European  bees  failed  to  become  established,  and  apparently  failed  to 
interbreed  with  the  indigenous  subspecies,  A.  m.  scute/fata  (Fletcher  1978). 
This  conclusion  may  be  supported  by  the  detection  of  only  one  variant  common 
(D404)  to  the  USA/European  and  South  African  samples.  The  existence  of  this 
variant  may  represent  a  relic  of  a  previously  successful  introduction  of  European 
genes  into  the  South  African  honey  bee  gene  pool,  or  it  may  represent  an 
ancestral  form.  ^         ^ 

A  more  exhaustive  survey,  including  the  identification  of  variants  and 
alleles  present  in  the  original  temperate,  European  populations,  may  indicate 
how  well  the  samples  of  the  present  study  represent  the  variation  throughout 
the  parental  and  neotropical  populations.  While  most  of  the  variants  and  alleles 
detected  in  the  neotropical  bees  were  the  same  or  similar  in  African  bees,  there 
is  evidence  that  European  markers  have  been  retained  by  neotropical  bees. 

In  the  region  of  Brazil  into  which  the  African  bees  were  released  there 
may  have  been  a  greater  concentration  of  west  European  compared  to  east 
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European  bees.   The  relatively  constant  frequency  for  the  Msp\  variant  MSOl 

in  neotropical  bees  suggests  that  it  is  a  neutral  marker,  incorporated  as  a  result 

of  west  European  introgression  into  the  African  population  as  it  was  becoming 

established  and  has  been  carried  along  as  the  bees  have  migrated  north  (Lobo 

&  Kreiger  1992;  Smith  1991). 

The  release  and  spread  of  African  bees  has  been  disruptive  to  the 

practice  of  beekeeping  in  the  neotropics  (Cantwell  1974;  Goncalves,  Stort  & 

De  Jong  1991;  Michener  1975;  Roubik  1980,  1989;  Spivak,  Fletcher  &  Breed 

1 991 ;  Taylor  1 977;  Winston,  Taylor  &  Otis  1 983).   Migratory  beekeeping  and 

commercial  queen-rearing  operations  in  the  USA  will  be  most  severely  impacted 

if  African  bees  are  not  excluded  (Danka,  Rinderer  &  Collins  1987;  McDowell 

1984).    Identification  of  honey  bees  will  be  essential  for  the  segregation  of 

African  bees  from  the  USA  beekeeping  industry,  an  integral  component  of  the 

intense  agricultural  system.     The  markers  reported  here  will  provide  more 

thorough  and  assurable  information  about  the  incidence  of  hybridization  in  the 

neotropics,  and  will  provide  a  means  for  following  the  progress  of  the  African 

bees  in  the  USA  (Hall  1990;  Hall  &  Muralidharan  1989;  Smith,  Taylor  &  Brown 

1989). 
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